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Abstract 

We develop a new method of constructing an asymptotic series in powers of N~ 1 / 2 as 
N — > oo for the function of N arguments which is a solution to the Cauchy problem for the 
equation of a special type. Many-particle Schrodinger, Wigner and Liouville equations for 
a system of a large number of particles are of this type, when the external potential is of 
order O(l), while the coefficient of the particle interaction potential is 1/N; the potentials 
can be arbitrary smooth bounded functions. We apply this method to equations for N- 
particle states corresponding to the iV-th tensor power of an abstract Hamiltonian algebra 
of observables. In particular, we show for the case of multiparticle Schrodinger-like equations 
that the property of iV-particle wave function to be approximately equal at large N to the 
product of one-particle wave functions does not conserve under time evolution, while the 
same property for the correlation functions of the finite order is known to conserve (such 
hypothesis being the quantum analog of the chaos conservation hypothesis put forward by 
M.Kac in 1956 was proved by the analysis of the BBGKY-like hierarchy of equations). In 
order to find a leading asymptotics for the iV-particle wave function, one should use not 
only the solution to the well-known Hartree equation being derivable from the BBGKY 
approach but also the solution to another (Riccati-type) equation presented in this paper. 
We also consider another interesting case when one adds to the iV-particle system under 
consideration one more particle interacting with the system with the coefficient of the 
interaction potential of order O(l). It happens that in this case one should investigate not 
a single Hartree-like equation but a set of such equations, and the chaos will not conserve 
even for the correlation functions. 







1 Introduction 



The chaos conservation hypothesis is well-known in statistical physics. This hypothe- 
sis put forward by M.Kac [1] in 1956 for the case of classical systems has the following 
analog in the quantum case of the system of N bose-particles moving in //-dimensional 
space. Consider the /c-particle correlation functions [2] 

K t ktN (x 1 ,...,x k ;y 1 ,...,y k ) = 

J dx k+1 ...dx N ^ t N (x 1 , x k , Zfc+i, zjv)\f^(yi, y k , x k+1 , ...,x N ) (1) 

corresponding to iV-particle wave functions ^^(xi, ...,x N ) which specify states of the 
system and satisfy the iV-particle Schrodinger equation (xi, ...,xn £ R" are particle 
coordinates, t G R is the time variable). The quantum analog of the chaos hypothesis 
is the following. Suppose that at the initial instant of time t — the correlator (1) 
factorizes as iV — » oo, k = const as follows: 

H\ tN {x u x k ; yi, ...,y k ) -> ip t (x l )...ip t (x k )ip t *(y l )...ip t *(y k ), (2) 

where </?* is one-particle wave function such that / ofo;|<^*(a;)| 2 = 1. Then the property 
(2) holds for arbitrary time t as well. 

The discussed hypothesis can be justified for the case of the external potential 
of order 0(1) and the particle interaction potential of order 0(1/N). The mathe- 
matical proof has been obtained in [3]. The method of justification of the property 
(2) was based on the ideas of [2]: the iV-particle equation for the density matrix 
VPjy-^i, x N )^f%(y 1: ...,yjv) was integrated over last N — k variables and the chain 
of equations for lZ k N was obtained in a way analogous to the method of derivation 
of the BBGKY hierarchy found almost simultaneously by Bogoliubov, Born, Green, 
Kirkwood and Yvon for the classical case. It was also found in [2,3] that the function 

obeys the Hartree equation being of widely use in physics for studying quantum 
systems with a large number of particles. 

The property (2) has the following physical meaning in terms of mean values of the 
observables being operators acting in Hilbert space L 2 (R uN ). Consider observables 
An with kernels of the special form 

Po l 

A N (x 1 ,..,x N ]y 1 ,...,y N ) = ^—— £ 

p=l ly P- l<i 1 ^...^i p <N 

A ip) (x il ,...,x ip ;y il ,...,y ip ) 5(x i -y i ), (3) 

where Po £ N,A^ are kernels of operators acting in L 2 (R^ p ).The property (2) implies 
that the mean values of the observables A N in the state have the following limit 

as N — > oo: 

(V^An&n) = j dx 1 ...dx N ^(x 1 ,...,x N )(A N ^ t N )(x 1 ,...,x N ) -> 
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Po l , 

E- / ( P t (x 1 )...(p t (x k )(p tm (y 1 )...(p tm (y k )x 
p=i P' 

xA {p) (y 1 , ...,y p ;x u ...,x p )dx 1 ...dx p dy 1 ...dy p . 

In this paper we consider a new formulation of a problem. Namely, we discuss if the 
property (2) is valid when k also tends to infinity, for example, k = N. We also study 
whether the chaos property allows us to find limits as N — > oo of mean values of 
observables of a more general form than (3). 

To solve these problems, we construct an asymptotic formula as iV — > oo not for 
the correlators (1) but for the full iV-particle wave function obeying the initial condi- 
tion of a product of one-particle wave functions (remind that a number of arguments 
of the wave function also tends to infinity as the parameter of the asymptotic expan- 
sion tends to zero). We will see that the asymptotics will not factorize into a product 
of one-particle wave functions. Therefore, the chaos hypothesis fails if k — > oo in 
(1). Thus, when one makes an attempt to find mean values of general observables 
uniformly bounded with respect to N, one can't use the property of factorizing of the 
wave function, contrary to the case of the observables of the type (3). 

Consider the multiparticle Schrodinger equation 

d 

ih—^ t N {x 1 , ...,x N ) = 



N 



£ [-2Tn A ' + Uix ' ) + N 

j=l \ LUl / ly l<i<j<N 



V t N {x 1 ,...,x N ), 



(4) 



where U is the Planck constant, m is the particle mass, Aj = d 2 /dxf is the Laplace 
operator, U is the external potential, j^V is the particle interaction potential being of 
order 1/N. The asymptotic formula for the iV-particle wave function is then expressed 
not only through the solution to the well-known Hartree equation 



h 2 

A + WHx) 



<p\x) 



where W l is a self-consistent potential 

W\x) = U{x)+ Jv(x,y)\ V \y)\ 2 dy. 
One should also use the following system: 



dt v 1 



h 2 

-—A + WHx) 
2m V J 



uHx) 



+<p t (x) I dyV(x,y)(i P t *(y)u t (y)+v t (y)if t (y)), 



(5) 



(6) 
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n 2 

-—A + W\x) 
2m 



(7) 

v\x) 



+tp t *(x) I dyV(x,y)^ t *(y)u t (y) + v t (y) V t (y)) 

The system (7) can be formally obtained by the following procedure. One can write 
the system consisting of the Hartree equation (5) and the equation conjugated to 
it, consider the variation system for it and substitute the variations of ip and ip* by 
u and v that should not be conjugated. This variation system with independent 
variations of ip and ip* coincides with (7). The asymptotic formula for the iV-particle 
wave function is expressed through the solution to eq.(5) and through the operator 
transforming the initial condition for the Cauchy problem for system (7) into the 
solution to this Cauchy problem. 

We can consider equations of a more general form than the multiparticle 
Schrodinger equation (4). Namely, we can study equations for the functions ^jv G 
L 2 (X), where X is an arbitrary measure space. Such equations are of the form 

ij t % = NA N ¥ N , 

where operator An has a kernel of the type (3). Notice that the Schrodinger equation 
is a partial case of this equation. Asymptotic solutions to it that obey more general 
initial conditions than a product of one-particle wave functions are constructed in sec- 
tion 5. These asymptotic formulas implying the results on the chaos non-conservation 
for eq.(4) are proved in section 6. In section 7 we evaluate the corrections to the 
asymptotic formula. 

We can consider not only Schrodinger-like equations but also sets of such equa- 
tions. The generalization of the method which is applicable to such case is to be 
developed in sections 8,9. The technique of these sections can be also used when one 
considers the iV-particle systems interacting with an additional particle, so that the 
evolution equation has the form: 

d ( h 2 \ N 

ih-¥ N (y, Xl ,...,x N ) = [N ^-—A y +U(y)j +^V(^,y) 

N ( h 2 \ l 

+ E -^T- A * + U{xi) + - Vix^xj^iy^,, ...,x N ), 

i=l \ zm / ly Ki<j<N 



(8) 



The first term of the right-hand side of this equation corresponds to the motion 
of the additional particle of the mass M/N in the external potential U(y) (y is the 
coordinate of the additional particle). The second term corresponds to the interaction 
of the additional particle with the iV-particle system. If the first and the second terms 
were of orders 0(1) and 0(1/N) correspondingly, one could apply the technique to 
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be considered in sections 5-7 without modification. But the coefficients are N and 1, 
so another approach is needed. It happens that the Hartree equation (5) should be 
modified, there will be no longer a single Hartree-like equation, there will be a set 
of such equations, and the analog of the chaos hypothesis will then fail even for the 
correlation functions. 

We can consider the problem of chaos conservation not only for iV-particle wave 
functions obeying the multiparticle Schrodinger equation but also for other cases. 
Namely, one can investigate iV-particle density metrices obeying the iV-particle 
Wigner equation [13] or iV-particle density functions (probability distributions) which 
obey the multiparticle Liouville equation [11]. One can find [11,13] asymptotic for- 
mulas for these cases. These asymptotics are also products of one particle densities 
at the initial time moment, while there is no such factorization at time moment t. In 
section 10 we will consider the equations for iV-particle states corresponding to the 
N-th tensor power of an abstract Hamiltonian algebra of observables. We generalize 
a notion of a half-density (discussed in [11,13] for different cases) by introducing the 
notions of a half-density representation of an abstract Hamiltonian algebra and of an 
abstract half-density. There are following examples of the latter notion: 

(a) the square root of the iV-particle probability distribution [11]; 

(b) the square root of the iV-particle density matrix [13]; 

(c) the iV-particle wave function. 

We show, that when our asymptotic method is applied to the iV-particle half- 
density equation, the average values of general bounded observables is unambigu- 
ously defined by our approximation for the half-density. The results to be obtained 
in section 10 imply, in particular, the asymptotic formulas found for the cases of 
Schrodinger, Liouville, Wigner equations. 

2 Violation of the chaos hypothesis for the N-particle 
wave function 

l.In the previous section we have seen that if the property (2) is satisfied, one 
can replace the wave function ^ N by the product of one-particle wave functions 
(p t (xi)...(p t (x N ) in order to find a limit as N — > oo of the mean values of the observ- 
ables of the special form (3). Consider the problem if such replacement is valid for 
finding mean values of general observables uniformly bounded with respect to N, 
\\An\\ < C. What is necessary and sufficient condition for this replacement? The 
following lemma tells us that such condition is 

J dx 1 ...dx N \^ t N (x 1 , ...,x N ) - c t N if t (x 1 )...<f t (x N )\ 2 ^n^oc (9) 

for some number c l N € C, \c l N \ —^n^^ 1. 

Lemma 1. Let <& N and ^ n be such elements of L 2 (W N ) that (& N: &n) = 
(Vn,Vn) = 1- 
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1. Let An be operators acting in L 2 (R, uN ) and uniformly bounded with respect to 
N,\ \A N \ | < C. Let for some set of numbers c^ G C, |cjv| — >jv-><x> 1 

(10) 

Then the property A N § N ) -^n~>oo A implies that N , A N ^/ N ) ^n^oo A 

2. Let 

(11) 

for arbitrary set of operators An uniformly bounded with respect to N. Then the 
property (10) is satisfied for some number cn G C, |cjv] — »jv-><x> 1 
Proof. 

1. Denote X N = — ctv^tv- We have 

(V N ,A N V N ) - ($ N , A N Q N ) 

= c N (X N , A N <Z> N ) + c* N (<5> N , A N X N ) + (X N , A N X N ) -^jv^oo 

because ||Aw|| < C and ||A^v|| — >n->oo 0. 

2. Consider the sequence X N of the form 

X N = V N - $ N ($ N ,ty N ). 
Notice that (<I>jv,-Xjv) = and choose the following set of operators A N : 

A N w = X N (X N , w),we L 2 (W N ) 
that are uniformly bounded with respect to N: namely, 

\\A N \\ = \\X N \\ 2 <{\\* N \\ + \\<$> N \\) 2 = 4, 
It follows from eq.(ll) that 

(V N ,A N V N ) - ($ N ,A N $ N ) = \(X N ,V N )\ 2 -\(X N ,Q N )\ 2 = 

= \(X N , X N )\ 2 — >jV-»<X> 0, 

so that the property (10) is satisfied for cjy = ($n,^n)- As 1 — \cn\ 2 = (^n — 
cn&n, — cn&n) —*n^oo 0, one has |cjy| — >jv-><x> 1- Lemma 1 is proved. 

2.Let U, V be smooth functions U : -> R, V : x R v -> R bounded with all 
their derivatives, V(x, y) = V(y, x). Denote by W™(R V ) the following set: 

^(R") = {/ : R" - C|/ e WftR"), k = TToo} 

Consider the initial condition y? : W — > C for the Hartree equation (5) such that 

^° G ^ 2 °°(R"),y |vAa;)| 2 da; = 1. (12) 
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The following lemma is proved in [4,5]. 

Lemma 2. There exists a unique solution ip l G W^^R?) to the Cauchy problem 
for eq.(5). 

As the Schrodinger equation is the partial case of the Hartree equation, lemma 2 
implies the following corollary. 

Corollary. There exists a unique solution ^/ f N G W^°(H I/N ) to eq.(4) which sat- 
isfies the initial condition 

f° N (x 1: ...,x N ) = p ( Xl )...p> (x N ). (13) 
It occurs that the property (9) is not valid. 

Theorem 1. Let V(x,y) ^ for any x,y G TV. Then there is no such interval 
[ti,t 2 ] that the property (9) is satisfied for t G [ti, t 2 ] f or some number c l N G C. 

This theorem is a corollary of the more general statement to be proved in section 

6. 

3. Let us consider a heuristic method to derive the result of theorem 1. Suppose 
that for some function c f N the initial condition (13) evolve into the wave function 

c t N ip t (x 1 )...ifi t (x N ) + z l N (xi, x N ), (14) 

where | \zn\ \l 2 ^n^oo 0. One can then expect that the property of chaos conservation 
for the full iV-particle wave function is also valid when the initial condition for the 
Hartree equation (5) is shifted by the quantity of order 1/y/N: 

ip° ^p° + 8<p°/y/N. 

In order to retain the property of the norm of the wave function to be of order 0(1), 
let us choose the variation Sip to be orthogonal to ip°. As the Hartree equation 
contains not only but also </?**, the function (p l transforms as follows: 

+ + + 0(1/N), (15) 



where A 1 and B l are some linear operators acting in L 2 (R^). Applying the transfor- 
mation (15) to formula (14) and multiplying Sip by e ta , one obtains that the following 
iV-particle wave function 

c^On) + -^=(A t Sip°(x 1 )e ia + B t (Sip°y(x 1 )e- ia ) + ...) 
viv 



x ...(p\x N ) + -^=(A t Sp°(x N )e M + B^Spyix^e-™) + ...) + z^O&i, *n) (16) 
V W 



is also an asymptotic solution to eq.(4). At the initial time moment formula (16) 
does not contain negative powers of e ta , while at time moment t such powers arises. 
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For example, the contribution of the power e m is zero at initial time moment and is 
equal to 

a t N <p t (x 1 )...(p t (x N ) 

1 N 

+c t ;^=E^(^i)-^ t (^-i)5 t (^ o )*(x i )^(x m )..^(x iV ) (17) 

VN i=1 

at time moment t for some numbers a.%, c l N G C (the term with a.% arises from 
the coefficient that may be a-dependent). Because of the linearity of eq.(4), the 
iV-particle wave function (17) is also expected to be asymptotic solution to eq.(4). 
Notice that the function B t (Sip )* can be decomposed into two parts: one of them 
being proportional to </?* and another being orthogonal to </?*: 

B\5^y = &y + B t '(5i P °y 

The contribution of the term to eq.(17) can be involved to the first term of eq.(17), 
while another term can't be treated in this way. As the norm of its contribution is of 
order 0(1), we are faced with the difficulty: the iV-particle wave function being equal 
to zero at the initial time moment evolves into non-zero wave function. The only 
possible way to resolve the difficulty is to adopt the violation of the chaos hypothesis 
(9) when B f ^ 0. The theorem 1 is heuristically justified. 

One can also expect that investigation of the operators A 1 , B 1 being obtainable 
from the variation system (7) can lead us to the correct asymptotic formula for the 
wave function. This is to be done in the following sections. 



3 Multiparticle canonical operator and asymptotic 
formula for the iV-particle wave function as N — ■> 

oo 

1. In the previous section we have seen that the wave function (17) may play an 
important role in constructing the iV-particle wave function asymptotics as N — > oo. 
We can also notice that such wave function satisfies the chaos property (2) and does 
not satisfy the propery (9). Therefore, the wave function (17) gives us an example of 
the state which can be replaced by the product <f t (xi)...if t (xN) in order to find limits 
of mean values of observables of the special form (3) but not of the general form. 

Let us give a generalization of the example (17). Introduce a notion of a multipar- 
ticle canonical operator being a partial case of the canonical operator corresponding 
to Lagrangian manifold with complex germ in Fock space [6,7]. 

Let us introduce the following notations. By T we denote the space of sets 
(go, gi(x\), #2(^1, X2), •••) of functions g n : H un — > C which are symmetric with re- 
spect to Xi e R^, belong to the space L 2 (R^ n ) and satisfy the condition 

00 . 

/ dx x ...dx n \g n [x x , ■■■,x n )\ 2 < 00. (18) 

n=0 J 
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By g n we denote the n-th component of g G T . Define an inner product in T as 

oo „ 

(g (1) ,g {2) ) = E J dxi-dx n g^*(x 1 ,...,x n )g^(xi,...,x n ). 



n=0' 



By ip G L 2 (R^), we denote the subspace of T which consists of all the elements 
g G T such that 

J dx 1 ...dx n (p*(x 1 )g n (x 1 , ...,x n ) = 0, n = 1, oo. (19) 
Consider the following element of L 2 (Tt uN ): 

N V 7 ? 

{K%9)(xu-,xn) = Y,-7= 9 P {xhi-,Xi P ) 11^)' ( 20 ) 

p=0 ViV l<i!<...<ip<JV j^l; 

where g G JF^. 

Remark. The wave function (20) satisfies the chaos property (2). 

Definition l.An operator : T v — > L 2 (R^ 7V ) o/ £/ie /orm (20) will be referred 
to as a multiparticle canonical operator. 

The following lemma shows that the norm of the function (20) is of order 0(1), 
although there are N\/(p\(N —p)\) terms in the sum over ii, i p , while the coefficient 
of this sum is of order N~ p / 2 , not of order 0(N~ P ). 

Lemma 3. The following relation 

N AH r 

\\K"g\\ 2 = N p( N _ p )\ J dx 1 ...dx p \g p (x 1 ,...,x p )\ 2 (21) 



is satisfied. 

Proof. One has 



p=0 ViV^ i< il< . ,.<i q <N l<j!<...<j p <N 

J dx 1 ...dx N g*(x il ,...,x iq )g p (x h ,...,x jp )Y[<p*(x i ) ]J <p(xj). (22) 

It follows from eq.(19) that for all non- vanishing terms in the sum (21) p = q, while sets 
ii, ...,i p and ji, ...,j q coincide. As the number of sets ii, ...,i p , 1 < i\ < ... < i p < N 
is equal to N\/(p\(N —p)\), eq.(21) is satisfied. Lemma 3 is proved. 

Remark. We can notice that the main contribution to the inner product 
(Ky g, K£ g) is given by the terms of the sum over p in eq.(20) which numbers are of 
order 0(1), not of order O(N). 

Corollary 1. | \K%g\ \ 2 < (g, g). 

Corollary 2. \\K^g\\ 2 -^jv->oc (g,g). 



Corollary 3. Let /, g be such non-zero elements of T v that for some cat G C 

\\K N J -c N K*g\\^ N ^ (yo 0. 

Then for some c G C / = eg. 

2. It occurs that the approximate solution to eq.(4) which satisfies the initial 
condition (13) should be expressed not only through the solution ip* to the Hartree 
equation but also through the solution to another equation. This is a Riccati-type 
equation: 

in-R\x, y) = V(x, y)^{x)^{y) + ( -—A x - —A y + W\x) + W\y) j R\x, y) 



+ 



J dy' V \y)V(y, y') V u (y')R\x, y) + J dx ip\x)V(x, x') V u (x)R\y, x) 

+ J dxdyR t {x ) x')R t {y,y')V{x\y)^{x)^{y) ) (23) 

where W l has the form (6). 

By W v we denote the space of complex functions R(x, y) : R" x R" — > C such 
that 

(i) ReWZ°(R 2 »),R(x,y) = R(y,x); 
(ii) 

J R{x,y)<p*{y)dy = -<p{x); (24) 
(iii) the operator M in L 2 (W) with the kernel 

M(x,y) =R(x,y) + ip(x)ip(y) (25) 

satisfy the property ||M|| < 1. 

Lemma 4. Let R° G W^o. Then there exists a solution R l G W v t to the Cauchy 
problem for eq.(23) with the initial condition R° . 

In order to prove this lemma, we will express the function R l through the evolution 
operator transforming the initial condition for variation system (7) into the solution 
to this system. 

Let us first construct this operator and study its properties. Denote by L, Y l and 
Z l the following operators in L 2 (R U ): 

L = - A A , (Z<w)(x) = k' 1 J dyV(x : y)<p t (x)<p t (y)w(y), 

(26) 

(Y*w)(x) = k~ 1 [W t (x)w(x) + J dyV(x,y)if t (x)ip t *(y)w(y)}. 



9 



The solution to the Cauchy problem for the system (7) is then as follows: 



u t \_ A 1 B l \ u° 
v f I l B u A u v° 



(27) 



where A\ B t are operators in L 2 (W) that satisfy the following equation 



.d_ ( A 1 B l 
% Jt\ B l * A 1 * 



L \ ( Y* Z l 
-L ) + 1 -Z** -Y** 



A 1 B l 
B t* A t* 



(28) 



and the initial condition 

A = E,B° = 0. 

We are going to show that the solution to eq.(28) is the following: 

where operators A* n , B l n are determined from the following recursive relations: 

( K Bt \ ( E \ 
\B» A£ ) [ E ) 

(A^ B^ n \ . /"* , / Y? ZJ \ / A T n _Y B^_i \ (or\\ 
bi: a* )-~ l Jo dT { -zr -rr ) { bzi ) ' 1 ' 

Notice that the series (29) is well-defined, since the functions U(x) and V(x,y), as 
well as the operators Y*, Z\ are bounded. The estimation 

• * || c n t n ,, , ,, c n t n 
K <— p Bl <— - 

for some constant c can be justified by induction. Thus, the series (29) converges. 
The following property is satisfied for operators (29). 

Lemma 5. The kernels of operators B t and A t R° belong to the space W£°(TL 2l/ ). 
Proof. Introduce notations: 

||5||= sup ||S£||,||S|| 2 = VfVS^S, 

H5ll=i 

||5-||M = ||(-A + l) m S(-A + l)- m ||, 
\\S\\t ] = ||(-A + l) m ,S(-A + l) m || 2 
for each operator S in L 2 (R U ). Let us show by induction that 

lKll (m) <^7-,ll^||5 m) <^7- ( 31 ) 
10 



for some constants c m . Inequalities (31) are correct if n — 0. Suppose them to be 
correct as n < k and check eq.(31) as n = k. It follows from eq.(30) that 

\\Al\\™ < f dr[\\Y;\\^\\AU\\^ + \\ZT\\<T ) \\BJ t _ 1 \t ) ], 

J 

< r^tlir/ll^ll^ill^ + II^II^II^Lill^]. 

JO 



t\\(m) 



\ B k\\2 



For c m > max[||17||( m >, ||^||^ m) ], eq.(31) is justified. 
It follows from eq.(31) that 



EiKii (m) <oo,Eii^ir<oo. 

n=0 n=0 

Therefore, 

WB'W^ < oo,|L4'||( m ) < oo, (32) 
WA'R ^ < WA'W^WR ^ < oo, 

Lemma 5 is proved. 

Let us show now that the operator (29) is really the solution to eq.(28). 
Lemma 6. The following relations are satisfied; 

( At+St _ At \ 

" i - ■ - LA* - (Y t A t + Z'B'*) )R% 



5t 



LB 1 — (Y t B t + Z t A t *)\\ 2 -^st^o 



11 St 

Proof. It follows from eqs.(29),(30) that 

At+5 [- At + iLA t + %{Y t A t + 2* B ») = ^ - 1 + iLSt A t 
St St 

+ (1 - e-^iij 1 A 1 + Z'B 1 *) 

_ e -iL { t +St) / dr^Y^A^ + Z^ B \r* - Y}Al - Z\B»). 
St Jo 

Making use of eq.(32) and of the relation 

e~ lLT -1 = [ T dae- iLa {-iL), 
Jo 

we find that the first statement of lemma 6 is satisfied. The proof of the second 
statement is analogous. Lemma 6 is proved. 

Let us construct the solution to eq.(23). 

Lemma 7. 



11 



1. The operator A 1 * + B l *R Q is boundedly invertable. 

2. The kernel of the operator 



R l = (B f + A'R^iA 1 * + B t *R°)~ 1 

belongs to W^t. 

3. The following relation is satisfied: 

T3t+5t pi 

||i n {Z 1 + (L + Y t )R t + R\Y l * + L)+ tfZ*#)\\ 2 (33) 

ot 

Proof. 

1. To prove the first statement, notice that eq.(30) implies that the matrix (29) 
is boundedly invertable: 

( # & Y 1 = A -n Bt n \( e^ \ 

where A l _ n , B l _ n obey the following recursive relations: 

( A Y Bt n \_ i f dr ( A T _ n+1 BL n+1 \ ( Y7 ZJ \ 



-n —n 



Bt. At„ I Jo 1 BL* n+1 XL* n+1 ; I -z? -rr 



which imply that 

( B u \ j ~ \ -Z? H Y r 



A t B t \ 1 / A t+ _ B tT 



(34) 



This means that the matrix (29) is a matrix of a canonical transformation [8]. This 
implies [8] that the operator A 1 is boundedly invertable. As 

E - ((A t )- 1 S*)((.4*)- 1 S t )* = (Ay^A 1 )*- 1 > 0, 

the operator (A'*) _1 _B** has the norm lesser than 1. Therefore, the operator 

(A u + B^R )- 1 = {E + (A t *)~ 1 B t *R°)~ 1 (A t *)~ 1 

is bounded, since = 1 and ||(A**) _1 5**|| < 1. The first statement of lemma 7 is 

proved. 

2. Let us check the property R* e W£°(R 2u ). It follows from the Cauchy-Schwarz- 
Bunyakovskii inequality that 



\\A M R t A K \\ 2 2 = Tr(A K R t *A 2M R t A K ) < 1 1 A 2M i?'| | 2 | | A 2K R l \ 
The first statement of lemma 7 and lemma 5 imply that 

\\A 2K R% < HA^^ + A^llalK^ + S**^)- 1 !! < oo. 
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Therefore, R* G W^(R 2u ). 

Let us check that i?'y?'* = —</?'. As the set of functions 

u* = i<p\ v* = -iip u (35) 

is a solution to eq.(7), one has 

(p** = - 5'V° = (A 1 * + B^R )^ *. 

Therefore, 

Eq. (24) is proved. 

Let us prove that the operator with the kernel (25) has the norm lesser than 1. It 
is sufficient to prove that (f , (E - R^R^O > 0, (£, (E - R t+ R f )0 = if and only if 
£ = Ay?'*. We have 

(£, (£ - R^R'K) = {{A u + B u R Q y 1 ^ (E - R 0+ R°)(A t * + B^R )' 1 ^). 

This quantity is non- negative and equals to zero if and only if (A t * J r B t *R°)~ 1 £ > = Ay? *, 
i.e. £ = Ay?'*. The second statement is proved. 

3. The proof of statement 3 is by straightforward substitution. Lemma 7 is proved. 

Remark. As the quantity h(Z f + (L + Y^R 1 + R t (Y t * + L) + RtZ^R?) coincides 
with the right-hand side of eq.(23), lemma 4 is a corollary of lemma 7. 

3. Let us give an asymptotic formula for the iV-particle wave function being a 
solution to eq.(4). Let y? G W^°(R V ),R G W^. Denote by the following element 



of Ty. 



®R,2n{xi,...,X 2n ) = 7=== E M (x il ,X ia )...M(x i2n _ l ,X i2n ), (36) 

2 n n\J(2n)\ i<n^...^ 2 „<2n 



®R,2n+l = 0,71 = 0,1, 2,3,... . 

where M has the form (25). The propery (19) is satisfied for $ R because of eq.(24). 
It can be shown by making use of the second quantization technique [8] that the 
property ||M|| < 1 implies that the series (18) converges. Therefore, G T 9 . 
Consider the solution to eq.(4) that satisfies the initial condition 



*° N = K^RO 



[N/2] 1 



E TT^Kin E M (x il ,Xi 2 )...M°(xi 2l _ 1 ,Xi 2l ) J] <P°(xi), (37) 

1=0 ) l<h^...^i 2 i<N i^i s 

where 

M°(x,y) = R (x,y) + cp (x)v (y),R G W>. 



13 



Consider also the solution tp* to eq.(5) which is equal to tp° at initial time moment 
and the solution R l to the Cauchy problem for eq. (23). Consider the functions: 



14/ 



(38) 



where 



H (<p*, <P) = jrJ dx<p*(x) (-^A + U(x) S j V (x) + 
^JdxdyV(x : y)\<p(x)\ 2 \<p(y)\ 2 ; 
c* = exp (—^f dr J dxdyV(x,y)ip T *(x)ip T *(y)R T (x,y) 



(39) 



Theorem 2. The following formula is satisfied; 



0. 



This theorem is a corollary of the more general statement to be proved in section 
6. It follows from theorem 2 and lemma 1 that one can use the iV-particle wave 
function 

W2] 1 



c^ NS \K^ Rt )(x 1 ,...,x N ) = c t e mst £ 



to 



x M t {x il ,x i2 )...M t {x i2l _ 1 ,x i2l )H<f t {x t ) 

l<h^...^i 2l <N i+i 3 

(M t (x,y) = R t (x,y) + tp t (x)tp t (y)) 



(40) 



instead of the exact wave function ^ N in order to find limits as iV — > oo of mean 
values of the observables uniformly bounded with respect to N. We see that not the 
form of the product of one-particle wave functions but the more complicated form 
(40) of the iV-particle wave function conserves under time evolution. The product 
(13) is a partial case of the wave function (40) which is realized when 



&(x,y) = -<p t (x) ( p t (y). 



(41) 



Therefore, one can make use of theorem 2 for finding an approximate solution to the 
Cauchy problem for eq.(4) with the initial condition (13). As the solutioin to eq.(23) 
does not, in general, have the form (41), the asymptotic solution to this problem 
has the form (40) with M t ^ 0. This implies that eq.(9) is not satisfied because of 
corollary 3 from lemma 3. Thus, theorem 1 is a corollary of theorem 2. 
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4 Heuristic derivation of the asymptotic formula 



In this section we consider a heuristic method to derive eq.(40). This method is 
analogous to the procedure of section 2 which is based on shifting the solution to 
the Hartree equation and allows us to conclude that the product of one-particle wave 
functions is not an asymptotic solution to eq.(4). 

As the chaos property (2) for correlation functions is satisfied for the iV-particle 
wave function (20) for arbitrary g G T^, as well as for the function ty f N [3], it is 
reasonable to look for the asymptotic expression for in the following form: 

V N = e iNSt K»g\ & 6^6 L 2 (IT), <7< G J> (42) 

It occurs that the heuristic method to be developed allows us to find the function 
S t up to an additive quantity that does not depend on the solution to the Hartree 
equation, while vector g l G T^t is defined up to a multiplier c* G C, \c*\ = 1. 

1. The function S* can be found by the following technique. Consider the small 
shift of the function ip by the quantity x/N of order 1/N in eq.(20). The function x 
can be decomposed into two parts: 

x = <p(<p,x) + x, (43) 

where x is orthogonal to (p. It follows from eq.(20) that the contribution of shifting 
by x to formula (20) is small. We will see that the coefficient of x chould be of order 
l/y/N for making such shifting of (p by x appreciable. The contribution of the first 
term of eq.(43) to eq.(20) is as follows: the p-ih term of eq.(20) is multiplied by 

(i + ^(<P,x)) N - p . 

As p = 0(1), all these quantities are approximately equal to exp(</?, Because of 
the remark after lemma 3, norms of terms of order p = O(N) in eq.(20) are small. 
Therefore, the following approximate formula takes place: 

K^ +j _ x g * e^K»g (44) 

This formula can be also proved rigorously [9]. 

Let us make use of eq.(44) for finding S f . Let us shift the initial condition for 
eq.(5) v? — > y?° + so that the solution to eq.(5) will be transformed as 

The factor S* depending on the initial condition tp° will be changed as S 1 ^ ) — > 
S\>p + ±x°)- It follows from eq.(44) that 

e iJVS*( V > +^X ) e (lP«,X«)-(v .X°) ~ e iNS*( V °)_ 
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Therefore, the variations of S* and y?°, </?* are related as follows: 

8S t = i[(<p t ,5 V t)-(<p°,8<p )]- (45) 

Eq.(45) determines the function S* up to a constant that depends on t but does not 
depend on y? . Let us check that eq.(38) really satisfies eq.(45). One has: 



55* = f dt f dx[i(8( f t )*(x)% t (x) - t-(8 V t )*(x)5 V t (x) - (8^)*(x 

JO J (Mb (Mb 



Scp* t (x) 



8(p t (x) 

The first term in this formula vanishes because of the Hartree equation, while the 
second term is equal to eq.(45). Thus, the function S* (38) can be found by the 
developed technique up to a function depending only on t but not on the solution to 
eq-(5). 

2. Let us derive conditions on the function g l e J 7 ^. We are going to consider the 
variations of </?* which are of order l/\fN and to obtain formula analogous to eq.(44) 
for the ./V-particle wave function 

K"iM + 0(l/VN)), (46) 



where g G T^. Notice that the vectors g + 0(1/ vN) belong to different subspaces 
JF , j_ of the space T and, therefore, should be different in general. 

First of all, consider the expression (46) in the case of vacuum vector g = 
i.e. go — 1, g a — 0, a > 1. The wave function (46) is then equal to the product of 
one-particle wave functions 

Mai) + 4=xM)~(<pM + -y^x(x N )) (47) 



which has been considered in section 2. Examine the case (ip, x) = 0. It follows from 
eq.(20) that formula (47) can be written as 

where C x is the following element of T^. 

C x , n {xx,...,x N ) = -i=xM-x{x n ). (48) 
yn\ 

When (ip, x) 7^ 0, one has 

K + ^ m = (1 + ^ X)fK«C x ,, (49) 
16 



where x is expressed from eq.(43). 

Consider now the following case of vector g in formula (46): g + 0(l/y/~N) = 
C fN J N = f + 0(l/VN).As 

(V + ^=X, fs) = 0, (50) 

one has 

K ^x C f» = K + ^( X+ M^ 0) = (1 + X + M) N K^ +fN -^ fN y (51) 

It follows from eq.(50) that 



(<P, /at) = -(X, f N )/VN = 0(1/ VN), 

(i + -j=(v, x + fN)) N *(i + x))^-^). 



Since (</?, /) = 0, one has 

K + ^ x {C f + 0(l/ViV)) = constK^e-^'^C^ + 0(1/VN)), (52) 

where constant in the right-hand side of eq.(52) does not depend on /. As the vectors 
(48) make up a full system of the vectors in the space [8], eq.(52) can be used for 
constructing an approximation for eq.(46). 

It is convenient to introduce creation and annihilation operators [8] in Fock space 
T . These operators are the following: 

1 ™ 

(a + (x)g) n (x 1 , ...,x n ) = —=^5(x - Xi)^ n _i(xi, ...,Xi-i,x i+1 , ...,x n ), 

(a~(x)g) n - 1 (x 1 , = y/ng n (x,X!, ...,x n -i). 

Notice that the operators 

a + [\] = J a + (x)\(x)dx,a~[\] = J a~(x)X*(x)dx, 

transform the subspace into if and only if (ip, A) = 0. Consider the operator 

V ViX = exp( J dx(a+(x)x(x) - a~(x)x*(x)), (53) 

where 



a [p (x) = a (x) — ip(x)a [ip], a*(x) = a + (x) — ip*(x)a + [ip]. 
The canonical commutation relations 

[a^x), a ± (y)] = 0, [a~{x), a + (y)] = 5(x - y) 
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and the representation of C/ through the creation operators 

C f = exp(a+[/])$(°) 
imply that the operator (53) transforms the vector (48) as follows: 

V VtX C f = exp(-( x ',/) - \( X ',x))C f+x >. (54) 
The factor e~ (x ' x >' 2 can be involved into the constant in eq.(52). Therefore, one 

has 

K^ x (9 + 0(1/ y/N)) = constK^ (V v , x g + 0(1/ VN)) (55) 

for the case g = Cf. Since the system of vectors (48) is full [8], eq.(55) is correct for 
arbitrary g. 

3. Let us find vector g t entering to eq.(42). Make use of eq.(55) for finding operator 
W* transforming g° into g t ,g t = W t g°. Suppose that variation of this operator is also 
of orde N~ l l 2 as the solution to eq.(5) is shifted by the quantity 0(1 /y/N). The 
iV-particle wave functions 

K^x* Wt 9 + 0(1/VN) 

and 

const(K^W%o iX og + 0(1/ y/N)) 

are then asymptotic solutions to eq.(4) if is a solution to eq.(5) and the set u f = 
= x l * is a solution to the system (7). As the asymptotic solutions coincide at 
initial time moment because of eq.(55), the same property should be satisfied at time 
moment t. Making use of eq.(55), one finds that 

V^tW* = constW%o :X o. (56) 

It is eq. (56) that allows us to find operator W* up to a multiplicative factor. It follows 
from eq.(27) that eq.(56) is satisfied if and only if 

J dx[a+ (x)(A t u° + B t v°)(x) - a - t (x)(B t *u° + A t *v°)(x)]W t = 

W J dx[a+ (x)u°(x) - a- (x)v°(x)\. (57) 

Eq.(57) is a straightforward corollary of eq.(56) when u° = x°i v ° = X°*- I 11 order to 
prove eq.(56) at arbitrary u°,v°, one should write eq.(56) for u° = x°G ta ,v° = -y^*e~ %a 
and integrate it over a with the weight e ±ia . 

Eq.(57) shows us that the operator W* corresponds to a linear canonical transfor- 
mation of the creation and annihilation operators [8]. It has been shown in [8] that 
W* is defined up to a multiplicative factor. In order to find the vector W t( 3>_R0 and to 
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show a role of the Riccati-type equation (23), it is convenient to introduce a notion 
of complex germ analogous to [10,6]. 

4. Let R G W v . Consider the following subspace Q R of the space L 2 (R U ) x L 2 (R U ) : 

Qr = {(v,u)\u = Rv}. 

Definition 2. A subspace Qr will be referred to as a complex germ corresponding 
to Re W v . 

The vector <&r satisfies the following interesting property. 
Lemma 8. 1. Let g G be such vector that 

J dx[a+{x)u{x) - a~(x)v(x)]g = (58) 

for any (v, u) G Qr. Then g = c$ R . 

2. The vector g = § R G satisfies eq.(58). 

To prove this lemma, one can make use of the definition of creation and annihi- 
lation operators, rewrite eq.(58) in terms of components g n of g and find them by 
induction to be equal to expression (36) up to a multiplier. 

Lemma 7 implies the following statement. 

Lemma 9. 

Q Rt = {(B'V + i'V.iV + BVjj^V ) G Q R o} 

Corollary. 

W*$ fl o = c% R t (59) 

for some constant c* G C. 

Namely, consider the vector W'fijo. It follows from eq.(57) that this vector satisfies 
the condition (58) for (t>*,-u*) G QRt, and, therefore, is equal to c$ R t. 

Therefore, theorem 2 is heuristically justified. Notice that the developed method 
can be applied also to the simpler case of the complex germ approximation for quan- 
tum mechanics as h — > 0, see appendix A for more details. 

Notice also that eq.(57) allows us to construct another asymptotic solutions to 
eq.(4) with the help of the complex germ creation operators: 

Kv = J dx[a;(x)v*(x) - a-(x)u*(x)}. (60) 

It follows from eq.(57) that 

W'A£ r .A^.4 R . = ^Zv A Zi*«- («) 

As any element of the space can be presented as a linear (maybe, infinite) combi- 
nation of vectors (61), eq.(61) allows us to reconstruct the operator W* and to find 
asymptotic solutions 
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to eq.(4). 

The approach developed here can be used for constructing asymptotic solutions 
to equations of a more general form than (4). Let us formulate the corresponding 
theorem. 



5 General case : formulation of the theorem 

Let X be a measure space, be a self-adjoint operator in L 2 (X N ) of the form 

H N = NHq + Hi + ... + N 1 ~ k H^ r , (62) 

where 

(H l N ^)(x 1 ,...,x N ) = Y / d yh--- d yi p Ki{x il) ...,x iv ]y hl ... ) y ip ) 

p=l iV V- \<i l ^...-Li v <N J 

x fy(x 1 , ...,^-1,^,^+1, ...,x ip _i,y ip ,x ip+ i, ...,xjv), (63) 

where H, (p) = # z (p) is a kernel of the operator acting in L 2 (X P ) which is symmetric 
separately over x^, ...,x ip and over y ix , ...,yi p and 

h\ p) *(xi, x p ; yi, y p ) = Hj p \y 1 , y p ; xi, 

Analogously to section 3, denote by JF the space of sets (go, gi(xi), 52(^1, ^2), •••) of 
functions g„ : A" n — > C which are symmetric with respect to Xi G X, belong to the 
spaces L 2 (X n ) and satisfy eq.(18). By C T we denote such subspace of T that 
consists of all the elements of T which satisfy eq.(19). The multiparticle canonical 
operator of the form (20) is denoted by : — > L 2 (X N ). Denote by $r e 
eq.(36). 

Let be a solution to the Cauchy problem 

ij^N = H N & N , 

(64) 

n=^A^...A^ flP , 

where A + has the form (60). 
Denote 

H/(y?*, V 5 ) = Y — / dxi...dx p dyi...dy p Hl P> (xi,...,x p ;yi,...,y p ) 
P =i P- J 

x ^*(xi)...^*(x p )^(yi)...^(y p ), G L 2 (X), (65) 



20 



Hi(ip*, ip) = H; (ip*, ip). Let </?* be a solution to the Cauchy problem for the following 
equation: 



such that / <ix|v9 (a;)| 2 = 1, i?* satisfy the equation: 

+ /»-rf(x,,') M ^ (y) + /d,Wrf(».»') M g^ (y) tf(l>',»). (67) 

where arguments V 7 * °f the function iJ are omitted, satisfy the following 

system: 

% Jt U ^ x) = J dy W(x)5<p(y) Ui{y) + 5^xW¥) VM )' 

By 5* we denote expression (38), while 

/ „t 1 /■ 5 2 H 

exp ( - Vo rfr 2 y ^ mx)m/ (i ' y) + Fi 



(68) 



(69) 



where arguments </?**, </?* of the functions H , Hi are also omitted. 
Let the following functions of xi, ■■■,x m , z±, z s - 2 j-f 

J dx m+1 ...dx p dy 1 ...dy p ip t *(x m+1 )...(p t *(x p )Hj; p \x 1 , ...,x p ;j/i, ...,%,) 

X(p\y s+1 )...(p\y p )R t (y 1 ,y 2 )...R t (y2 j -l,y2j)R t (y2j+uZ 1 )...R t (y s _ h Z s _2j-i) 

( 7 °) 

belong to L 2 ^™-*- 2 *'-*), where Ji, J* e 1, fc, X j = v} - Ru*j. 
Theorem 3. T/ie following relation takes place: 

||^ _ cV^T^A^V^+^ll -^iv^oo 0. 

Remarks. 

1. For some special choice of H&\ eq.(4) is a partial case of eq.(64), the Hartree 
equation (5) is a partial case of eq.(66), eq.(23) is the analog of eq.(67), eq.(68) is a 
variation system (7). Therefore, theorems 1 and 2 are corollaries of theorem 3. 

2. An asymptotic formula being approximately equal to the wave function 
accurate to 0(N~ L l 2 ) for arbitrary L > is to be presented in section 7. 
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3. The requirement for the functions (70) to belong to L 2 means that the N- 
particle wave function (61) belongs to the domain of an operator H N . This require- 
ment can be easily checked for the case of bounded operators H{* , while for general 
case one should prove this independently. For the case of conditions of theorem 2, the 
square integrability of eq.(70) is a corollary of the properties of the solutions to the 
Hartree equation (5) and to Riccati equation (23); some other cases are presented in 
[11-13]. 

4. We have used different notations, and h\ p \ for the same quantity. This 
has been done in order to simplify formulas to appear in section 8. In that sectoion 
we will denote by an operator in Z 2 , while its eigenvalue will be denoted as H^. 
The same remark is correct for the notations and 



6 Proof of the theorem 

1. Let us consider the more convenient representation for the operator Hn that allows 
us to find a commutation rule between and multiparticle canonical operator . 
Consider the space L 2 {X N ) as a subspace T N of the Fock space T of the form 

T N = { g ef\g a = 0,a^N} 

Consider the operator 

H i* = S J dx 1 ...dx p dy 1 ...dy p U\ p> (x 1 ,...,x p ;y 1 ,...,y p ) 

x a + {xi)...a + (x p )a~(yi)...a~(y p ). (71) 

Lemma 10. The operator (71) transforms F N into F N and coincide on the 
subspace J rN with the operator (63). 

The proof is by making use of the definition of creation and annihilation operators. 

As the operator Hn is expressed through the operators a^, it is sufficient to find 
their commutation rules with the multiparticle canonical operator. 

Lemma 11. The following relations are satisfied: 



a^(x)a , rJ N _ N 



K=Ka+{x), 

(72) 



a-{x)K% = ^K%a-{x). 



The proof is by making use of the following expression for the element g e T N C 

N 



If 

K*g = J2^=r / dx 1 ...dx p g p (x 1 ,...,x p ) 
p=0 VP- J 
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(i„'„(xi)a [ip] a+(x p )a [ip] a H 



and of the commutation relations between operators a ± (x). 
Consider now the operator jj.K^ . 
Lemma 12. The following relation takes place; 

|^ = <( | +0+Ma; ^ ] + (w | p)(w _ ft) 

+ (74) 

where an index t on ip is omitted, 

°J[x] = a+ M - a + [(p]((p : x):a~[x] = a~[x\ - a' [<p](x, <p), n = J dxa+(x)a~(x). 
Proof. Consider the iV-particle wave function 



d d " HvT 



dt ^ dt NPpl (N - p)\ 
x / dxi...dx p g p (xi, ...,x p )a + (xi)...a + (x p )a + [(p] N ~ p &°\ g G T^. (75) 



When one takes a derivative with respect to t, there will be two terms: one of them 
contains ^g p , another contains J^a + [</?]. Consider the first term. The function ^ can 
be decomposed into two parts: 



dt^ p ^ X ^' "' 



The second part being equal to 



belongs to J 7 ^, since the property g G J 7 ^ (19) conserves under time evolution, and 

V dt' 



contribute to eq.(75) as K£ ^-g. The first part gives rise to the following contribution 



to the expression (75): 



N I TV 1 1 r _ d 

xa + {x 1 )...a + {x p _ 1 )dx l ...dx p - 1 (a + [ V ]) N - p+1 <S>^ 
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which can be presented as 



K»(N-h)a- 



dip 
~dt 



Consider now the terms containing -|a + [ip\. They contribute to eq.(75) as follows: 



p=0 



NPpl (N - p) 



J 9 P (xi, ■.■,Xp)a + (x 1 )...a + (x p )dx 1 ...dx p 



d 



x(iV-p)a + [-v,](a + M) Ar - p $ (0) . 



According to the definition of a^, this expression is equal to 

Combining all the terms, we obtain eq.(74). Lemma 12 is proved. 
Corollary. 



.d_ 
~dt 



- H N y» st K» = e^ty (i± + ia^]a- t [j/] - 7{ N ) , (76) 



where 



d „ t . t d t . . n . . /——; 1 r d 



d 



dt 

k 



dt' 



N' 
pi 



v L dt 



N ^ L df 



4+7 - J dx 1 ...dx p dy 1 ...dy p 



Po v 

h P =inio m\(p -m)\8\(p-8)\ 

xip t *(x m+1 )...ip t *(x p )n^ ) (x 1 , ...,x p ;y!, ...y p )ip t (y s+1 )...ip t (y p )a+ t (x 1 )...a+ t (x m ) 

x (1 - n/N)(l - (n + l)/iV)...(l - (n + p - m - 1)/A0<y (y s ). (77) 

Remark. The corollary implies that the operator H' N is a product of iV by a poly- 
nomial in N' 1 / 2 : 

H' N = NH' + N 1 ' 2 ^ +H' 2 + ... + N x ~ k I 2 H k (78) 
for some K. The coefficients i?Q, if^, can be presented as follows: 

(79) 



dx 



<5(/9*(a;) dt' 



'5(p(x) dt 



(80) 
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6<p*(x) 

~\ I tedwWM 5 J x ^ {y) + j dxdy[\a+ t (x) — ^ {y) a+ (y) 

+ fl > (l) ^S(/^ (!,) + ^^ (l) M*)My) a;t(y)1 ' (81) 

where the arguments y?**, tp* of the functional H are omitted. One can notice that 
eqs.(38),(66) provide the nullification of H' ,H' V 
2. Let us prove that 



i^-^e^^'H-^ooO, (82) 



d 
It 

where 

*-<?«&■■<«**■ 

Lemma 13. 

\\H' k g t \\<oo,k = 0J{. 

The proof is by making use of commutation relations between creation and anni- 
hilation operators, formula 

= exp Q J dxdya+{x)M{x, y)a+(y)) $< > (83) 

and the condition of the square integrability of the functions (70). 

As H = H 1 = 0, for checking eq.(82) it is sufficient to prove the following lemma. 
Lemma 14. 

(i- + ia + [p]a-[- ( p]-H' 2 )g t = 0. 
Proof. 1. It follows from straigthforward calculaton that 

[ij t + m"V] V ~ K A&+ ]/« = 0,fe (84) 



if and only if 



5(p*(x)5(f*(y) 
6 2 H 

dt v ' ' 7 w \^</?(a;)&/?(y) vy/ ' S(p(x)S(p*(y) 
where 

i/ = u * - y? V, u*), / =v t - 



25 



a*, (3* are some complex functions. This system for u* ,v f , as well as eq.(84), is 
satisfied when (u t ,v t ) satisfies eq.(68). 

2. It is sufficient then to prove lemma for the case k — 0. From eq.(83) one has: 



d d dc* i f dM t 

l?^** = + 1 y dxdy(M t (x,y) - vW(y)) ) 

+| / dxdya% t (x)aj (y) [^ ( ^ (y) + / ^ ( "' ^ 



Lemma 14 is then proved as a corollary of eqs.(67),(69). 

Eq.(82) implies the statement of the theorem,the proof is analogous to [14]. Con- 

•N 

d 



sider the quantity e lNS * K^g 1 — being equal to 



e iNS 



l K4 -** N = fj Te -^^\i± - H N ){e iNST K^ - W N ) 
The following estimation takes place 

iNSt K^ - nil < f Q d ^t ~ H N )e wsr K»gr\\ 0. 



e 



Theorem 3, as well as theorem 1,2, is proved. 

7 Corrections to the asymptotic formula 

Let us construct now the iV-particle wave function that approximates the solution to 
the Cauchy problem for eq.(64) accurate to 0(N~ L ^ 2 ) for arbitrary L > 0. It happens 
that such wave function has the form: 

¥ NjL = e^K^g* + N-Wg{ + ... + N'^gU), 

where vectors g l m G T^t are of the form 

din = J2 77 / dxx.-.dxng* (xi,...,x n )at(x 1 )...at (x n )® R t, 



no(m) is a finite quantity, functions g m , n {xi, x n ) are to be defined by induction. 
Let g % m n be defined for m < I and all functions 
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x i P t (y s +i)-^ t (y P )R t (yi, y 2 )-R t (y2j-i, y 2j )R t {y 2 j+i, z 1 )...I?(y a - n , z s - 2 j- n ) 

) (85) 
belong to £ 2 (,Y m - s - 2 .?-™) f or m < I. Define the functions Xi n ( x ii ■■■■> x n) being sym- 
metric with respect to Xi and obeying the requirement / dx^* {x 1 )x t l ^ n {x 1 , ...,x n ) = 
by the following relation: 

X! ~~j=\ J dxi..dx n xl n (x 1 ,...,x n )a+t(xi)..a+ t (x n )$i i t = 



= H' 3 gl 1 + ... + H K gl K , (86) 

where all = by definition as n < i. Note that definition (86) is correct and 
Xi^ n = for sufficiently large n, the sum in the left-hand side of eq.(86) is then finite. 
When I = 0, let xl,n = by definition. 

Let g\ n be a solution to the Cauchy problem for the following equation: 



d d\nc\ t " f r 5 2 H 

<dt-^r )g ^ Xl -> Xn) = fJ dyk 

k=l 



S(p*(x k )5(p(y k ) 



I 



ot f ^ H o i*/ (m + l)(m + 2) 



5(p{z k )5(p{y k ) 2 

x J dy 1 dy 2 j^-^^gl n+2 (yi,y2,x 1 ,...,x n ) + x[ n ( x i, -, x n), (87) 

where c* has the form (69). Let be a solution to eq.(64) that satisfies the initial 
condition 

= K°(9°o + N-^g° + ... + N'^-^gU). 
Theorem 4. Under the conditions of theorem 3 

\\^n-^n,l\\ = 0(N~ l /'). 

Remarks. 

1. The requirement for the functions (79) to be square integrable is analogous to 
the same requirement for the function (68) and means that the asympyotic solution 

l belongs to the domain of an operator H^. 

2. For the case of eq.(4), the solvability of eq.(81) and the square integrability of 
eq.(79) can be proved in a way analogous to the proof of solvability of eqs.(7) and 
(23). When H N is a bounded operator, such statements can be also proved. 

Proof. It is sufficient to show that 

Wi-Hs)& N A\ = 0(N- L l% 



dt 
i.e. 

(ij t + ™V\ %Aj t ^ - H '^ m = H'^ln-i + - + H K gl_ K . (88) 
The proof of this relation is straightforwrd. Theorem 4 is proved. 
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8 Some aspects of problems with operator- valued 
symbols 

We have seen that the discussed method of constructing asymptotic solutions can 
be applied to equations of the form (64). However, one can be interested in the 
problem of generalization of the considered approach to the case of the set of such 
equations. Investigation of it is very important when one considers the quantum 
mechanical system consisting of two subsystems (one of them is the examined system 
of N bose-particles, another subsystem interacts with the first one), some examples 
are to be discussed in section 9. Notice that some of the results to be obtained in this 
section and in section 9 can be also derived by making use of the technique analogous 
to the derivation of the Ehrenfest theorem in ordinary quantum mechanics [20], see 
appendix B for more details. Let us now consider the specification of the form of the 
set of equations which is to be approximately solved. 

Let X be a measure space. Denote by I 2 <g> L 2 (X N ) the Hilbert space of sets of 
complex functions ^i(xi, Xn), I = 1, oo, x±, G X such that \& i G L 2 (X N ) 
and dxi...dxN\^i(xi, ...,xn)\ 2 < oo. Let H N be a self-adjoint operator in I 2 <g> 

L 2 (X N ) of the form (62); operators have the form (63), where ^ el 2 ® L 2 (X N ), 
while n\ p \xi, ...,x p ;yi, ■■■,y p ) being operators in I 2 with matrices 
Hyij(xi, ...,x p ; yi, ■■■,y p ) (I, J — 1, oo) are kernels of the following operators in 
I 2 ® L 2 (X N ): 

(Hi p) ^) I (x 1 ,...,x p ) = J dy 1 ...dy p H$ J (x 1 , ...,x p ;yi, ...,%,)* j(yi, -,y P ), 

* G l 2 (g)L 2 (X N ), 

which are also required to be self-adjoint. Denote by I 2 <g> JF the space of sets 
g = (9o,i, 9i,i(xi), 92,i(xi, x 2 ), ■■■) of functions g nJ : X n -> C, I = 1, oo, which are 
symmetric with respect to Xj G X, belong to the spaces L 2 (X n ) and satisfy the 
following condition analogous to eq.(18): 

oo oo „ 

dx 1 ...dx n \g nJ {xi,...,x n )\ 2 <oo. 

1=1 n=0 J 

By I 2 ® C I 2 ® J 7 , where tp G L 2 (X), we denote, analogously to the case of section 
5, such subspace of I 2 <g> T that consists of all the elements oil 2 ® T that satisfy the 
condition 

J dx 1 (p*(x 1 )g n j(x 1 ,...,x n ) = 0. 

Denote by K^f : I 2 <g> — > I 2 ® L 2 (X N ) the following analog of the multiparticle 
canonical operator (20): 

(^S)i(xi,-,x N ) = (K"g I )(x 1 ,...,x N ), 
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where g G l 2 ®^, while gi is the element of T of the form (go,i, gi,i(xi),g2,i(xi,%2), •••) 
at fixed /. 

We are going to find approximate solutions to eq.(64) for the case of the operator- 
valued function and ^ N G I 2 ®L 2 (X) by the technique analogous to the method 
discussed in sections 5-7. These asymptotic solutions are to be looked for in the 
following form: 

^ Ar = e^ t Kj(g*+7V-V2 g t + ... )) (89) 

where g* G I 2 <8> J-^t. We will formulate the theorem in section 9, while in this section 
we are to find g* heuristically. One should substitute the expression (89) to eq.(64) 
and make use of the commutation rule between operators and (id/dt — H N ). 
Lemma 14. The following relation is satisfied: 

{ij t - H N )e iNSt K» = e* NSt K»(iD t - <h' N ), (90) 

where H' N has the form (77) and 

D t = d/dt + a + [ V t }a- t [j tV t }. 

The proof of this lemma for the case of the operator-valued function is anal- 
ogous to the proof of eq.(76) for the case of sections 5-7. 

By (ip*,ip) we denote the operator in I 2 of the form (65). The operator 7i N 
in I 2 ® is then written in a form (78), where operators H Q , H 1 , H 2 have the form 
(79), (80), (81). For the simplicity, the operators in I 2 of the form like XE, where A is 
a number, are denoted by A. 

One can notice that one should choose g^, g*, ... in such a way that 

(_jV#; _ N l l 2 H[ + iD t -H' 2 - N^ 2 H' 3 - ...)(g* + N- l ' 2 g\ + ...)= 0{N~ L ' 2 ) (91) 

when the asymptotics accurate to 0(N~ L ^ 2 ) is looked for. 

An interesting feature of the operator-valued case is that the operators H , H 1 

cannot be set to zero by varying </?*, since the operator H (<^*,v 5 ) is not, in general, 

equal to H ((p*,(p). Therefore, in order to provide satisfaction of the relations like 

eq.(82), one cannot choose g* to be independent on N, g* must be choosen as g* = 
g * +A r-i/2 g t +iV -i g t. 

Let us consider the recursive relations for g*, which are derivable from eq.(91). 
First of all, consider the term of order O(N) in eq.(91) which has the form i? go = 0. 
It follows from eq.(79) that H' can be presented as H (y?**, y?') — A* for some number 
A*. Therefore, g(, should be chosen as gg = ( ® gg, i.e. 

9o,n,l( X l> X n) — (l9o,n( X l> •••> X n), 

where ( G I 2 is the eigenvector of the operator H (y?**, acting in I 2 . 
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Let Holf 1 *, ip l ) be eigenvalue of the operator Ho(y2'*,y?') and smooth function of 
t, n(<^'*,<^ t ) be projector on the corresponding eigenspace, so that 

(H (yA ^ - H (^, ^))I%<*, = 0. (92) 

The term of order O(N) vanishes then, if S* has the form (38). 

Suppose that this eigenspace is one-dimensional. The method under consideration 
can be also applied in analogous way to the case of finite and (^'-independent dimen- 
sionality of the eigenspace. The case of terms intersection, ehen this dimensionality 
depends on ip, requires the more careful treatment. 

Assume that H is an isolated point of the spectrum of the operator H in I 2 , so 
that there exists a unique operator R such that 

iffl = 0, (H - H )R(1 - IT) = 1 - n. 

We will denote this operator R as 

(H - Ho)-\l - n) = R, 

the arguments tp**, <p l of the operator II and functional Hq are omitted. 

To each operator A in I 2 with the matrix Ajj we assign the operator in I 2 <g> JF^ 
that transforms the vector with components g n j(xi, ...,x n ) into the vector 
Yfj=\Ai,jg nt j(xi, ...,x n ). This operator in I 2 <g> T v will be also denoted by the same 
symbol, A. 

Let us consider other terms of eq.(91). It is convenient to present g* as 

Si Si Si i 

where 

g! ll =iig|,gf L = (i-n)g|. 

It follows from eq.(38) that H' gf = 0. Therefore, eq.(91) can be written as 

(H - H )& + + (H' 2 - iD t )gl_ 2 + ... + H K gl_ K = 0. (93) 

The vector g^ is determined in a unique fashion from this relation if and only if 

IT^g^ + (H' 2 - iD t )gl_ 2 + ... + lt K f^ K ) = 0, (94) 

since n(H — Hq) = (because H is a self-adjoint operator). 

When m — 1, eq.(94) implies that HH^q = 0. It follows from eq.(80) that one 
should require that 
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i.e. 



/. d t( , <5H 

n i—ip (x) — - — — — 
V dr v ; 6<p*{x) t 



n = o. 



(95) 



This equation is the analog of Hartree equation for the operator-valued case. 

Furthermore, one can find g^-i from eq.(93), substitute it to eq.(94), make use 
of formula (95) implying that HH^g^^ = and obtain the equation for g t m _ 2 - In 
order to find its solution, one can first find g^~_ 2 from eq.(93) and obtain the following 
equation for g„_ 2 : 

-UH'^Ho - H )~\l - + g™- 2 ) + (H' 2 - A)^_ 3 + # 3 <4-3 + ...) 

+ U((H 2 - <A)(gJ!-2 + S™- 2 ) + •••) = 0, (96) 
where g„_ 2 are found from eq.(93). Therefore, one can look for the quantities g^~ 
and gjjj by induction. Let go" 1 , go", gm_ 3 , gm-3 t> e already found. Then one should 
find g„_ 2 from eq.(93) and g„_ 2 from eq.(96). The equation for g„_ 2 has the form 
like 

U(iD t -H' 2 + ^(Ho - H )-\l - n)i/;)gl 2 = xt-2 (97) 

for some right-hand side Xm-2 sucn that (1 — H)Xm~2 = 0- Eq.(97) is analogous then 
to eq.(88). 

Let us simplify our main equations, (95) and (97). As I1(H — H ) = (H — i^o)n = 
0, one has 







dx 



5 (n(H.-*)n) = n(i^)n. 

5(H - H ) 



5(f*(x) 

Therefore, eq.(95) takes the form of eq.(66), while 



a + (x) 



'ggo - Hp 

5(p*(x) 



+ a (x) 



S(p(x) 



In order to simplify eq. (97), denote by G I 2 , ||C*|| = 1, the eigenvector of H (y?**, </?*)■ 
One has 



g2!- 2 = C ®9 l m -2^x 

Making use of the relations like 



t\\ 

m-2 



Cm-2- 



(9* 



n 



*(H - gg) 

5<f(x) 

} 5U 6{H -H ) 

J S(p(y) 6<p(x) 



5U 



5cp(x) 



(H - if ), 



n + i/ (Ho ~ go) n = o, 

5ip(x)5ip(y) 



5ip{y) 5ip(x) 
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6U U) 6(H -H ) I[ _ 5U 6(H -H ) u 



5tp(y) 6<p(x) 6<p(y) 6<p(x) 

and commutation relations between operators a ± (x), one obtains that eq.(97) takes 
the following form: 

f 1 5^H 
[Wt + lt ~J ^ [ 2< t(x) ^(x)V(y) < <(y) 

+ < t(x) ^^(y) V (y) + ^^ (g) Mx)My) a ^ (y)lg "- 2 = x ™" 2 ' (99) 

where 7* is a number of the form: 

Notice that eq.(92) has been already solved in sections 6,7. Therefore, the functions 
g l m are found. 

Note that the leading asymptotics has the following form: 

V N = gV^KX ® a£+...a£+ <£>**) + 0(1/VN). (101) 

1 ' 1 k ' k 

It is remarkable that the functions wjj., -R* obeys the equations coinciding with 
the equations obtained in sections 5-7 for the case of a single Schrodinger-like equation, 
not for a set. The only difference with the considered case is that the quantity (100) 
arises in the left-hand side of eq.(99), so that the phase factor c* differs from the 
factor c* obtained in sections 5-7.This confirms the heuristic arguments of section 4 
that predict the form of the asymptotics by making use of the equation for tp* only. 

Consider the terms of the additional factor 7* in more details. The first term was 
studied in ref.[14] for the case of semiclassical approximation (the Maslov canonical 
operator with real phase) for quantum mechanics. When one considers adiabatic 
perturbation theory, this term known as Berry phase [18] also arises. 

Consider the last term of eq.(100). One can formally rewrite it as 

1 f J ,„ i 2 (H„ -#„),,, 



ij^ M,^.] ^ (io2) 

We will study some examples for which the form (100) is correct, while the expression 
(102) contains divergences to be eliminated. Therefore, we will use eq.(100). 

Discuss now the problem of chaos conservation for the operator-valued case. One 
can study different eigenvalues and eigenvectors of the operator H and obtain asymp- 
totic solutions being superpositions of the formulas like eq.(101): 
00 

E cW NSt K 7 ((j ® 7 -< + J , w £ **,) + o(i/Vn). (103) 

j=i ' ' j> j' 
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One can consider the quantities analogous to the /c-particle correlators (1) 

Tll^ix!, ...,x k ;yi, ...,yk) = 

oo . 

53 / dx k +i...dx N i& t N j(xi, ...,Xk,x k +i, ...,x N )V%j{yi, ...,y k ,x k+1 , ...,x N ), (104) 
j=i j 

where G I 2 <S> L 2 (X N ). The correlation functions (104), as well as the correlators 
(1), allow us to predict the limits as N — > oo of mean values of the observables of the 
special form (3). We can notice that for the element (103) of the space I 2 <8> L 2 (X N ) 
such correlators (104) have limits 

oo 

Xk',yi,-,yk) -»• TV^oo 

j=i 

where Aj are some numbers. For example, one can choose the functions tp\ to be 
coinciding at the initial time moment. Since the functions ip\ obey different Hartree- 
like equations corresponding to diferent eigenvalues of Ho, they will not, in general, 
coincide at time moment t. Therefore, the chaos property (2) being satisfied at t — 
does not hold at arbitrary time moment. Therefore, even for the correlation functions, 
the chaos conservation hypothesis fails for the operator-valued case. 



9 Operator- valued case: the theorem and some ex- 
amples 

1. Let H ((p*,(p) be non-degenerate eigenvalue of the operator H (y?*,y?), tp* be a 
solution to eq.(66), S* have a form (38). Let g n satisfy recursive relations (93) and 
the norm of the vectors H m g n , D t g n be finite. Consider the solution to the equation 

ij t ^ N = H N v* N , ¥ N el 2 ® L 2 {X N ) 

that satisfies the initial condition 

K = K£(gg + iV-^ g o + ... + AT-^gO). (105) 

Theorem 5. The following relation is satisfied: 

\\¥ N - e^K^g* + 7V-V2g* + ... + iV-^gDH = 0(N-^). 

Remarks. 

1. The initial condition (105) for the Cauchy problem is not arbitrary. In partic- 
ular, g° should be an eigenvector of the operator H . Moreover, the vectors g^~ can 
be expressed through g^_ l5 g° with the help of eq.(93). When the initial condition 



33 



does not satisfy these properties, one can present it as a superposition of the permiss- 
able initial conditions corresponding to different eigenvalues of H . Therefore, one 
should make use of the solutions to different Hartree-like equations (66). The chaos 
property (2) will not conserve then under time evolution. 

2. When the initial condition (105) is allowable, one can solve the Cauchy problem 
for the recursive relations by induction with the help of the technique analogous to 
the previous section: one can first express g^ through g^-i, ...,go, then one should 
find the solution to the Cauchy problem for eq.(97) by using the substitution (98) 
and reducing eq.(97) to eq.(99). One can notice that if the initial condition for g* 
is expressed as a result of action of a polynomial in creation operators o^(x) to the 
vector Q R , then the function Xm-2 * s a l so expressed in such a way: 



xl 

— '71! 



I^^r / dx 1 ...dx n xj, n (xi,...,x n )a\(x 1 )...a+ t (x n )$ R t. 
n Vnl J 

The vector function g^„ 2 has then the form 

gz = ex P(^ r * rft )H^y / dx 1 ...dx n gl n (x 1 ,...,x n )a+ t {x 1 )...a+ t (x n )<5> R t, 
where / = m — 2, 

r '^'-U d ^j£m M '^ < ioe > 

7* has the form (100), while the set of functions gj (xi, x n ) is a solution to the 
Cauchy problem for eq.(87). 

3. The proof of theorem 5 is analogous to the proofs of theorems 3,4. Nevertheless, 
one should take into account that the function 

K5(g*+ArVV 1 + ... + A^/ 2 gi). 



approximately satisfies eq.(64) accurate to 0(N~ L 2 1 ), not to 0(N^ L i 1 ). Therefore, 
one should substitute to eq.(64) the asymptotic formula with two additional terms. 

4. If the dimensionality of the eigenspace is a constant D more than 1, one should 
consider the orthonormal basis (I, ...,(b i n this eigenspace and present g„_ 2 as 

D 

gm-2 = Yl Cj ® Sm-2- 
1=1 

The equation for g^ 2 has then also the form (99), but 7* should be considered as a 
matrix D x D of the form 

t ( \f d C N \ ,m 1 f , , r / n ^ 2 Ho , x , J 2 (H -H ) 
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<fy>(ar) <5y?*(x) 



57V 



2. We have considered the case of finding asymptotic solutions to the infinite set 
of equations for functions ^i(xi, xn), I = l,oo. The same method is applicable 
when one considers the set of d equations for d functions ^i(xi, xjv), / = l,d,i.e 
one studies the equation for the element of R d <E> L 2 (X N ). 

Example 1. Let us consider the simple example. Consider the set of two 
Schrodinger-like equations 

i%— I ^n,i( x i,--,x n ) \ = y- h ( B u(xi) B l2 {xi) \ ( V t NA (x 1 ,...,x N ) \ 
1 dt { ^ t N o(x 1 ,...,x N ) ) £t I B 21 ( Xi ) B 22 {xi) ) I ^ t N Jx 1 ,...,x N ) ) 



+ 



N / *2 



i=l \ / iv l<j<j<iV 



^ N J Xl ,...,x N ) ) ^ 



This set of equations corresponds to the following physical problem. Besides N- 
particle system, there is a two-level system interacting with N particles. It is the first 
term in the right-hand side of eq.(107) that describes this interaction. One can note 
that the coefficient of this interaction is of order 0(1), contrary to the coefficient of 
the particle interaction potential which is 1/N. Therefore, the term with the matrix 
B is to give rise to the additional term in the Hartree equation. 

According to the developed technique, consider the matrix H (y?*, (p) correspond- 
ing to the operator in R 2 

2 / B n (x) B 12 {x) \ 
{ B 21 (x) B 22 (x) i ' 



H (</?*,</?) = H$(<p*,<p) ^ J l ) + / dx ^ x 

+1 kJ dxd y v ( x ^y)\^( x )\ 2 \^(y)\ 2 - 



The operator H has two eigenvalues: 

Hfa; <p) = # V, <p) + Hu + H22 ± /%*, <p), 

where 



H IJ (i P \ V ) = Jdx\<p(x)\ 2 B IJ (x),[3= ] J( 11 2 22 ) +H 12 H 21 . 
Therefore, there are two Hartree-like equations: 



&*v = M)- (108) 
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One finds the following asymptotics for the solution to eq.(107): 



t ^l>'-K^(C T f6\ { h *^) + 0(l/ViV), 

1,— ' 1,— fe_,— ' 

where S± has the form (38), T l ± has the form (106), (p± obey eq. (108), (ui t ±, i>i,±) 
obey the variation system (68), i?± obey eq.(67), 

s± u t± i^ 22 ~ i o\ 

d = aH 12 ,( 2 =a{ 2 ^' 

a =2[3{ ±13). 

Let us now consider the problem of divergences in eq. (102). When one formally 
calculates the quantity (C, s v *\x)slp{x) Q > one obtains: 

5(0)[^(S 11 (x)+ J B 22 (x)) 

±2^ ((S u (x) - B 22 (x)) Hu ~ H22 + tf 21 5 12 (z) + tf 12 5 21 (*))]. 

It happens that this infinite quantity is equal to the divergent part of (£, ^^jf^^ C) ■ 
Therefore, the divergences can be eliminated in eq.(102), but they arise in calculation. 
In order to avoid arising of infinite quantities, we have used eq.(100) instead of eq. 
(102). 

Example 2. The developed approach can be also applied to the more interesting 
case of eq.(8). To construct asymptotic solutions to eq.(8),one should first consider 
the operator 

Ho = H° + l (-^A y +U(y) + / dzV(z,y)Mz)\^ 

depending on ip and find the eigenvalues Hq\lp*, ip) and eigenfunctions (i[tp*, vMu) °f 
this operator. Then one should make use of the solutions to the Hartree-like equations 
(66), as well as to eqs.(67),(68). The asymptotic solution has then the form (103). 
The argumentation on the chaos non-conservation for the correlation functions which 
has been presented in the end of section 8 is also valid for this example. 
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10 The problem of chaos conservation and 

asymptotic formulas for abstract Hamiltonian 
algebras 

1. We have constructed asymptotic solutions to equations of the type (64) and to sets 
of such equations. It was shown that the L 2 (X N ) norm of the difference between exact 
and approximate solutions tended to zero as N — > oo. Since was interpretted as 
a multiparticle wave function, lemma 1 told us that its approximation e tNS K^tg 1 
constructed in sections 5,6 could be used instead of for finding limits as iV — > oo 
of mean values of general observables uniformly bounded with respect to N . 

On the other hand, there are some other interesting cases. For example, one can 
consider the case of classical statistical mechanics when the algebra of observables is 
presented not as an operator algebra but as algebra of real-valued functions on the 
phase space. The case of quantum statistical mechanics when states are specified not 
by wave functions but by density matrices can be also studied [13]. Examination of 
these cases requires one to formulate analogs of lemma I. 

In this section we are going to generalize lemma 1 for the case of abstract Hamilto- 
nian algebras of observables (see, for example, [3,19]) which involves all the cases con- 
sidered earlier. Our results will imply, for example, the conclusion of refs. [11,13] that 
asymptotic solutions to Liouville and Wigner equations which are found by the de- 
veloped technique should be interpretted not as approximate densities but as approx- 
imations for iV-particle half-densities which are equal to the square roots of density 
functions (matrices) and also obey Liouville (Wigner) equations. We will introduce 
the notion of an abstract half-density which generalizes the notions of refs. [11.13]. 
We will also consider the equation for it and find its asymptotic solutions. 

These asymptotics are to be expressed through the solutions to eqs.(66),(68). If 
we consider and </?** to be independent, the set of equation (66) and equation 
conjugated to it will form a Hamiltonian system playing an important role in con- 
structing tunnel asymptotics [13]. It happens that for the case under consideration 
in this section such Hamiltonian systems can be simplified, so that the number of 
equations can be cut in half. 

2. Definition 3. [3,19] A Hamiltonian algebra A is a set of complex linear space 
A and mappings n : A x A —> A, \ : A x A —> A, j : A —> A denoted also as 
n(A,B) = AB, \{A,B) = {A, B}, j(A) = A + , A,BeA, if the following axioms 
hold: 

Al ). for any A, B,C E A,a, (3 <E C 

a) A{aB + (3C) = aAB + 0AC, {A, aB + (3C} = a{A, B} + (3{A, C}; 

b) (A+)+ = A, (aA + (3B)+ = a*A + + (3*B + ; 

c) {AB) + = B+A+, {A, B} + = {A + , B + }; 

d) {A, BC} = {A, B}C + B{A, C}; 

e) A(BC) = (AB)C; 
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f) {A, B} = — {B, A}; {{A, B}, C} + {{B, C}, A} + {{C, A}, B} = 0; 
A2). there exists such element I G A that AI = I A = A, {A, 1} = for any 
A g A; 

A3), for any A,B G A and some % G R AB - BA = ih{A, B}. 
Remark. The case % = is also allowable, so that one cannot write {A, B} = 
UAB-BA). 

By C we denote the complex linear space of all linear functionals p : A — > C. 
Introduce also the notation 

C + = {p G £|p(J) = 1, P (A + ) = (p(A)r,p(A + A) > OVA g £}. 

Suppose that £ + ^ 0. 

By Ciwe denote the set of such elements A G A that sup pe£+ |p(A)| < oo. 
Consider the following functional p : Af — * R: 

p(A) = sup \p{A)\. 
P ec+ 

Lemma 15. The following relations are satisfied: 

p(A + B) <p(A)+p(B),p(aA) = \a\p(A), A, B G A F , a G C. 

The proof is straightforward. 

Let {A G Af\p(A) 7^ 0} 7^ 0. Denote by Cp the linear space of such linear 
functionals Af C that sup p( - yl ) =1 |p(A)| < oo. Inroduce the following norm in C F : 

HpII = sup \p(A)\, p G C F . 

p(A)=l 

Lemma 16. The functional \ \ ■ \ \ : C F — > R satisfies the following properties: 
l)||ap|| = |«|||p||, ||pi + p 2 || < IIPill + ||P2||,||p|| > 0,p,pi,p 2 e C F ,a e C; 

2)\\p\\ = 0^ p = 0. 

Proof. The proof of the first property is straightforward. Let us prove the second 
property. When p = 0, is obviously equal to 0. Let ||p|| = and show that p = 0. 
One has: p(A) = for any A G A F such that p(A) 7^ 0. It is sufficient to check that 
p(B) = when p(B) = 0. Lemma 15 implies that p(A + B) < p(A) + p(B) = 
p(A),p(A) < P(A + B) + p{B) = p(A + B). Therefore, p(A) = p(A + B) ^ and 
p(B) = p(A + B) — p(A) = 0. Lemma 16 is proved. 

Remark. The Hamiltonian algebra A playes the role of the algebra of observables. 
The elements of C + specify possible states of the system. The quantity p(A) is then 
the average value of the observable A in the state p. The observables of the form A + A 
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are usually called non-negative. The property p(A + A) > for elements of C + means 
that average values of non-negative observables are also non-negative. The functional 
p(A) has the following physical meaning: it is the largest possible average value of the 
observable A. The role of the norm is the following: if the quantity \\pi — yo 2 1 1 is 
small, the difference between average values pi(A) — P2(A) of any abservable A such 
that p(A) = 1 is also small. 

Let us give now some examples of Hamiltonian algebras. 

Example 1. Denote by A c n the algebra of smooth functions A(pi, q±, ...,p n , q n ) 
on H 2un . Define the algebra operations as follows: 

(AB)(X) = A(X)B(X),A + (X) = A>(X),X = {p 1 ,q 1 ,...,p n ,q n ); 

(109) 

, „ „ w x " fdAdB 8AdB\ . 

{A, B}( Pl , qi , ...,p n ,g n ) = g^_____j qi , ..., Pn , q n ). 

It is easy to see that the axiom Al is satisfied. The element I(X) = 1 playes the role 
of /, so that the axiom A2 is also checked. The axiom A3 is satisfied when U — 0. 
Therefore, the algebra A° n is Hamiltonian. Note that it is the algebra of observables 
for classical satistical mechancs. 

Lemma 17. The functional p(A) is the following in the case of example 1: p(A) = 
supxei&™\A(X)\. 

Proof. Let \A(X)\ < C. Then 2C - A(X)e^ - A+(X) e -^ = B+B for some 
observable B. Therefore, for any p e C + one has < p{CI) — Re[e 1 ^ p{A)) ,i.e. 
|p(A)| < C. Thus, p(A) < sup XeIt 2vn \ A(X)\. Consider now the element p G C + of 
the form p Xo {A) = A(X ). We see that p(A) > \p Xo {A)\ = \A(X )\. We prove lemma 
17. 

Let us give the examples of elements of C F . 

A) To any function p from L 1 (R 2un ) one can assign te element p G C F of the form 

p(A) = J dXp(X)A(X),X G R 2un . (110) 

One has: ||p|| = J dX\p(X)\. 

B) The following element of Cp'- px Q {A) = A(X ) can be also formally written 
in the form (110), but p(X) is a generalized function p(X) = S(X — X ). One has 
IHI = 1- 

Note that the function p playes the role of a probability distribution. 
Example 2. Let A q n be algebra of bounded operators in L 2 (R n ). Let 

ir(A,B) =AB,\(A,B) = ~z{AB — BA),j(A) = A + ,I = E. 

Analogously to the previous example, one can check axioms A1-A3. The functional 
p(A) is equal to the ordinary operator norm: 

p(A) = \ \A\ \ = sup (ip,Aip). 
IMI=i 
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Checking this property is analogous to the proof of lemma 17: one can present op- 
erators cE — Ae % ^ — A + e^ as B + B. Note also that elements of Cp play the role of 
density matrices. 

3. Let us now generalize the notion of half-density to the case of an abstract 
Hamiltonian algebra. 

Definition 4. A half- density representation of a Hamiltonian algebra A is a set 
of a Hilbert space Ti and mappings IT : Af x T~L — > Ti, A : Af x 7~t — > 7i, denoted also 
as 

U(A, <p) = A<p = Il A (p, A(A, <p) = {A, ip} = A A ip, A e A F , <p G W, 

i/ £/ie following axioms hold: 

HI), the mappings Yi A and A A are linear operators in 7i which are defined on a 
common domain T> C T~i; 

H2).U. 1 = E,A J = 0; 

H3). the following relations are satisfied for any A,B<E Af, <p, X G "P.' 
a J U aA+/3B = aTL 4 + (3U B , A aA+(3B = aA A + {3A B , 

b) {A,B<p} = {A,B}<p + B{A,<p}, 

c) U AB = ILTP 3 , A< A > B > = A A A B - A B A A , 

<*; (<p,u a x ) = (u A ^,x).(^A A x ) = -(A A +y?,x)- 

Elements ofTL are called half- densities. 

To any element of (f G V one can assign the following element of C F : 

Pv (A) = {<p,n A tp). (in) 

It is important that the norm \ \p ipi — p<p 2 \ \ is small when \ \ipi — tp 2 \\ is small. 
Lemma 18. The following relation is satisfied: 

\\p Vl ~ P<p 2 \\ < \\<Pi ~<P2\ 1(1 kill + II^H) (H2) 



Proof. One has 

Hpvi -Pvall = su p l(^i,n A v?i) - (^2,n A ^ 2 )| (113) 

p(A)=l 

Notice that p v G C + when \\<p\\ = 1. therefore, |p^(LI A ) | = | (c^, II" 4 ^) | < 1 in this 
case. Thus, ||vr A || < 1 when p(A) = 1. Eq. (113) implies then eq.(112). Lemma 18 is 
proved. 

Let p G Cp. Denote by {H, p} the following functional 

{H,p}(A)=p({A,H}). 

Definition 5. Let H G A,H = H + . An abstract Liouville equation is the 
following equation for p l G Cp, t G R: 

% = {H,P t }- (H4) 
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An abstract half- density equation for tp 1 ETC is 



% = {H^}- (H5) 

Lemma 19. Let ip* obey eq.(115). Then p v t obey eq.(114)- 
Proof. Let A G A. One has 

± P(pt (A) = (vMrW) + (AV,nV) = 
= [n A , a><) = (yAn^V) = /v({A #})■ 

Therefore, eq.(114) is satisfied. Lemma 19 is proved. 
Let us give examples of half-density representations. 

Example 1. Consider the Hamiltonian algebra A° n . Define oprations AB = 
U A B, {A, B} = A A B, where A G A c n , B G L 2 (R 2m ) by eq. (109). Let V = S(R 2un) . 
The axioms of the half-density representation are satisfied. We can also notice that 
the element p v G Cp is equal to 

p^(X) = \p(X)\ 2 . 

for example, if tp is real, p is presented as a square of the function p. Because of this 
reason, the function tp has been called half-density function in ref.fll]. 

Example 2. Consider the Hamiltonian algebra A^. Let Tt be Hilbert space 
L 2 (X n ). Determine operators U A ,A A as 

U A tp = A<p, L A p> = ^-Atp. (116) 
in 

It is not hard to check axioms Hl-H3.The half-densities play the role of the wave 
functions in this example. The density matrix p v is proportional to the projection 
operator on one-dimensional subspace. 

Example 3. Consider the same Hamiltonian algebra A q n . Choose the Hilbert 
space Tt as the space of Hilbert-Schmidt operators on L 2 (X n ) with the inner product 
{tp, x) — Trtp + x- Let n A , A A be the following: 

U A tp = Atp,L A tp = ^[A,tp]. (117) 

The density matrix P(p is Ptp = tptp + ; the matrix p has been called half-density [13] 
for reasons analogous to example 1. 

We can notice that different equations (Schrodinger and Wigner) are treated from 
the same point of view: they are abstract half-density equations for the cases of 
examples 2 and 3 correspondingly. Contrary to the half-density representations, the 
Hamiltonian algebras are identical for examples 2 and 3. 
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4. For the simplicity, consider the case when elements A of the Hamiltonian 
algebra A are presented as functions A Y on a measure space y, while the algebra 
operations are as follows: 

(AB) X = J dYdZd XYZ A Y B z ,{A,B} x = J dYdZf XYZ A Y B z , 

(118) 

(A + ) x = A* x , X,Y,Z ey,A,B e A, 

f and d are (generalized) functions y x y x y — > C. To simplify the notations, denote 
the integrals like (118) as dxYzA Y B z ,fxyzA Y B z , i.e. we integrate over repeated 
indices. 

Elements of Cf are presented as (generalized) functions on y, the quantity p(A) 
can be formally written as p(A) = pxAx- 

A tensor product .A®^ = A x ... x A is then presented as an algebra of functions 
A Yi ...y n : y x ... x y = y N — > C, Yi, Yjv G 3^. The algebra operations are defined 
as: 

AT 

{A, J B}x 1 ...x JV = XI dxiyiZi-"dx p _iy p _iZ p _i/x p Y p z p 
P =i 

X C?x p+ 1 Z p+ 1 Y p+ 1 • • • ^Xjv Zjv Yv ^Yi . . . Y N B Zl . . . Z N , 

(AB) Xl ...x N = d Xl Y 1 z 1 —dx N Y N z N A Yl ...Y N Bz 1 ...z N , (A + ) Xl ...x N = A* Xi _ Xn , 

A,Be A® N , Xi, Yi, Zi e y. 

It is not hard to check axioms of a Hamiltonian algebra for A® N . Note that the axiom 
A3 is important for such check. 

Let 7 be a polynomial functional 7 : C — > R of the form 

for some observables G A® k . 

Definition 6. [3] An e-uniformization of the functional (119) is a set of observables 
H (N) e A m . 

K Jfc-1 

k=l h - l<ji^...^j k <N p^ji 

Consider now tensor products of half-density representations. Without loss of 
generality, one can assume that the Hilbert space 7i is presented as L 2 (X) for some 
measure space X . Elements of Ti will be denoted as ip x , x G X; integrals like / dx(p*x x 
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will be denoted as (p^Xx- Let (L 2 (X),U, A) be a half-density representation of the 
Hamiltonian algebra A. The operations II, A are 

(Aip) x = bxxyAx^y, {A, (fi} x = C X xyA X Py 

for some (generalized) functions b, c:yxXxX^C. Consider the Hilbert space 
n ®N = Hx xH = L 2( X N} and the following operations: 

{Ap) Xl ...x N = bx 1 x 1 y 1 ---bx N x N y N Ax 1 ...X N fy 1 ...y N , 

(121) 

N 

xi. ..xn / j ^X\x\y\---bx v -\x v -\y v -\C-X v x v y v $X v +\x v +\y v +\---$XNXNyN 
p=l 

xA Xl ...X N P yi ...y N , 

where fixxy = bx xy — ifocx X y,Xi G y,Xi,yi G X. One can check the axioms of 
definition 4 for the operations (121). This half-density representation is called a 
tensor product of N half-density representations (L 2 (X), II, A). 

Let A" be fixed, e = 1/N, H have the form (120). Then the abstract half-density 
equaton (115) takes the form 



d ^° 1 k 

u,/ fc=i jy i<ii<...<« fc <A r g=i 



X Ps p+1 x lp+1 y lp+1 --Ps k x lk y lk fxi...y ll ...y lk ...x N - ( 122 ) 

Lemma 19 tells us that the element of Cp of the form 

Pxt..X N = {P<p*)xi_...X N = bx 1 x 1 y l ---bx N x N y N ¥x* 1 ...x N ¥ t y 1 ...y N 

obeys the abstract Liouville equation (114). 
Remarks. 

1. The algebras A C N and A q N considered above are the following tensor powers, 
A C N = (^l)®^; A C N = (AD® N . The N-th tensor powers of half- density representations 
(116) and (117) of the algebra A\ are half-density representations of the algebra A q N 
which are given by the same equations, (116), (117). The half-density representation 
of the algebra A C N which has been considered in example 1 is also the AT-th tensor 
power of the analogous representation of algebra A\. 

2. Eq.(114) has the form of the Liouville equation for the case of algebra A C N and 
of the Wigner equation for algebra A q N . For the case of the system of A" particles 
moving in the external potential U and interacting each other with the potential j^V, 
the observable H entering to eq.(114) has the form of uniformization (120) of some 
functional 7 which has the same form 



l[p\ = J dpdqp(p, q) + U(q) 
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- J dp 1 dp 2 dq 1 dq 2 V(q 1 , q 2 )p(pi, 9i)p(p2, <b) 



+ 2 

both for the classical and quantum cases. One should only take into account in 
quantum case that p(p,q) is a symbol of the operator p G C in L 2 (Tl u ) which is 
defined as 



1 r 



where px is a kernel of the operator p. 

3. For the case of the Hamiltonian (120), eq.(114) has been considered in [3]. It 
was shown that the property of the correlation functions 

T^x^.'xk = Pxf...x Ar / x fc+1 ---/x Ar (123) 
to be approximately equal to the products of one-particle correlators 

K^fi^N^ooP^-Px, (124) 
conserves under time evolution, while px obeys the abstract Vlasov equation 

This justifies the chaos conservation hypothesis. 

5. One can notice that eq.(122) for the half-density is of the type (64), if = 
NHq,Hq has the form (63) and 

K 



Ho(tf*,(p) — ^2 ^^^^■■■S k ^Sixiy 1 ---bs q - 1 x q - 1 y q -JCs qXq y q Ps q+ ix q+1 y q+1 --Ps N x N y N 



1 

k=l ^ 5=1 

x V*x 1 - ( P*x k l Pyi- l Pyk- ( 126 ) 
Therefore, the technique developed in sections 5-7 can be applied to eq.(122). One 
can construct such function (p^/ G L 2 (X N ) that 

Was -V iVas ~¥ ) -^N^oo 0, (127) 

where tp N,t is a solution to the Cauchy problem for eq.(122). Lemma 18 tells us 
that estimation (127) means that one can use approximate half-densities in order to 
find limits as iV — > oo of average values of the observables uniformly bounded with 
respect to N. In particular, one can confirm the chaos conservation hypothesis for 
the correlation functions of finite orders and deny it for the iV-particle densities in a 
way analogous to section 3. 

It is interesting that eq.(124) is the analog of the Ehrenfest theorem. Namely, one 
can introduce creation and annihilation operators a* 1 and apply lemma 10. Eq.(120) 
takes the following form: 

= H (a + /VN, a-/VN)<p»*, (128) 
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where tp* is such element of T that only A^-th component of it differs from zero, while 

K ° 1 k 

H (a+/VN, a'/VN) = £ ^ £ Wj*ls k 

k=l ' q=l 

x ^Six 1 y 1 ---^S q -ix q - 1 y q - 1 iCs qXq y q f3s q+ ix q+1 y q+1 --f3s N x N y N <lti--- a t k a y 1 --- a y k - (129) 

The correlation functions TZx^ k .!x k can be presented through the average values of the 
function of operators a ± / \J~N : 

~-,(N,k,t) (N — k)\ Nt , , _ - Nt\ 

K x i.: Xk = m {if ' ,b XlXiyi ...b XkXkyk at 1 ...at k a yi ...a ykV > ). 

As the commutator between operators a ± /\/ r N tends to zero as A" — > oo, one obtains 
eq.(125) from Heisenberg equations. Thus, one confirms the chaos conservation for 
correlation functions. 

Let us present the theorem for the A^-particle density which implies the results of 
[11,13] for multiparticle Liouville and Wigner equations. Let H have the form (120), 
e = 1/N. Consider the solution to eq.(114) that satisfies the initial condition 

PXi...X N = bx 1 x 1 y 1 ---bx N x N y N vi-vn 

Let (p f be a solution to the Cauchy problem for eq.(66), R l be such a solution to 
eq.(67) that obeys the initial condition R°. Denote 

& = exp -- / dr- — - — R T , 

V 2 Jo 5^ x 5^ y xy J 

Pxt a .X N = \ ct \ 2 bx 1 x 1 y 1 ---bx N x N y N (K^t^Rt)l 1 _ XN (K^Rt) yi ...y N 

Theorem 6. The following relation takes place: 

II N,t,as N,t\\ n 
WPn ~ PN II ~^N^oo u. 

This theorem is a corollary of theorem 3, lemma 18 and reality of S f (eq.(38)). 

Consider now eq.(66) for our case in more details. Note that it can be treated as 
a Hamiltonian system if one introduces the following Poisson brackets: 

{Vx, ¥y} = W* x , y?*} = o, {lfi x , y?*} = -%. 
One can present the functional H Q as 

H = ^\ rf \p]- 1 \p-ha]],h^0; 

(130) 
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where 

OX = ic Xx y<P* x <P*y, PX = b X xytp* x <P*y (131) 

The Poisson brackets between p, a are the following: 

{PX, Py} = hfzXYPZ, {PX, Oy} = fzXYPZ, W X , &y} = fzXY^Z- 

Therefore, the Hamiltonian system in terms of p, a is divided in quantum case into 
two independent parts: as {px,PY — Ti&y} = 0, the equations for p and p — ha are 
independent. In classical case the Hamiltnian system consists of two equations: the 
equation for p coincides with eq.(122), another equation for a is linear. 



11 Conclusions 

We have developed a new asymptotic method that allows us to find approximations 
for functions of a large number N of arguments as N — > oo, as well as the corrections 
to the leading order of the asymptotic formula. We have seen that this technique is 
applicable to eq. (64) being of a general form, as well as to the set (finite or infinite) of 
such equations. Multiparticle Schrodinger, Liouville and Wigner equations are partial 
cases of eq.(64). 

We have noticed that for the case of Schrodinger equation our approximate wave 
function can be used instead of the exact wave function for finding limits of mean val- 
ues of general observables uniformly bounded with respect to N. We have considered 
the case of a general Hamiltonian algebra, justified the chaos conservation hypothesis 
for the correlators and denied it for the iV-particle densities. It is interesting that for 
the operator-valued case the chaos does not conserve even for the correlators. 

Two methods have been used for constructing such asymptotics. One of them is 
heuristic and allows us to construct the asymptotic formula up to a multiplicative 
factor c t e iN , where constant C* does not depend on the solution to the Hartree-like 
equaion (66), while c* depend, in general, on this solution. We can see from theorem 
3 that the reason for this arbitrariness is as follows. When we change the operator 
iff, the Hartree-like equation (66) does not change, while the number c* is multiplied 
by a (^-dependent factor. The function C* will change when one adds the constant 
h l to the operator Hq . As the heuristic method of section 4 uses only the form of 
eq.(64), not the form of the hamiltonian, it cannot predict c*,C*. 

Note that the argumentation of section 4 is also applicable to the operator-valued 
case. This implies that the only difference of the asymptotic formulas for the ordinary 
and operator- valued cases is the value of the constant c*. This conclusion is justified: 
all quantities entering to eq.(lOl) are identical to the analogous quantities in the 
asymptotic formula of section 5, except for the phase factor c*. 
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The relation between the considered technique and complex germ method [10] 
can be investigated in more details [6,15-17]. As we have seen in section 10, one 
can consider such representation for the A-particle wave function that eq.(64) will 
transform into ordinary Schrodinger-like equation, while the analog of the Planck 
constant will be 1/N. The asymptotics constructed in this paper correspond [6,17] 
to the isotropic manifold [10] of the special form. General Lagrangian manifolds with 
complex germs can be also considered when one studies the infinite superposition of 
the obtained asymptotic formulas, see [6,15-17] for more details. 
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Appendix A 

As mentioned in conclusions, our asymptotic method for finding approximate solu- 
tions to eq.(64) as N — > oo is analogous to the complex germ technique. Namely, 
eq.(64) is analogous to the ordinary Schrodinger-like equation 



ih^^- = H(x,-ih-^]^ t (x),xeR n , (132) 



the Hartree-like equation (66) for complex function is an analog of the Hamiltonian 

system 

for two real vectors P*, Q l G R™. According to section 10, the chaos conservation 
hypothesis for correlation functions resembles the Ehrenfest theorem. The multi- 
particle canonical operator is analogous to the canonical operator for the complex 
germ in a point, the constructed asymptotic solutions resemble the wave-packet-like 
approximate solutions in quantum mechanics 

tf{x) = c t e ^ + pt {x - Q t)) g t (^S^j + 0(Vh), (134) 

g* E S(R n ). 

When initial conditions are given, the functions c*, can be found by solving 

the equations obtained in complex germ theory by subsitution of eq.(134) to eq.(132). 
The purpose of this appendix is to show how such equations can be heuristically 
derived without such substitution in a way analogous to section 4. 
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First of all, find the dependence of S* on P*, Q t . When one varies the initial 
condition for the classical trajectory by the quantity of order h, vectors P*, Q f varies 

as 

P* -> P* + MP*, Q* -> Q* + MQ*, 
the only change of the wave function (134) as H — > is multiplication by the quantity 

i(<SS*-P'<SQ«) 

where ^S** is a variation of the quantity S^P , Q°) depending on the initial condition 
for the classical trajectory. Therefore, 

5S l - P*5Q* = const, 

i.e. 

S t = f t [P T —Q T -H(P T ,Q T )]dT + C t (135) 
Jo dr 

for some quantity C* that does not depend on the classical trajectory. 

Consider the variation of the solution to eq.(133) by a quantity of order Vh: 

P t -> P* + v^P', Q* -> Q* + v^Q', 
The function (134) transforms then into the following function: 



c t e i ( ^ + P*(,-o*))^ fc-^ + o(Vh), (136) 



where 

/(0 = canst exp ^ ^P'f - ^J^)) 
As the functions (134), (136) should be asymptotic solutions to eq.(132), the operator 

SP^-SQ'--^- (137) 

should transform solutions to the equation for g l to solutions, when (SP*, 5Q V ) is a 
solution to the variation system 

d d 2 H d 2 H 

dt SQ = dPdP 6P + 8PdQ 6Q ' 



(138) 



~dt 6P = dQdP 6P + dQdQ 6Q ' 
Notice that (<5P*, 5Q*) may be complex quantities. 
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Therefore, if (5P°£ - SQ°\-^)g° = then 



(6I*Z - SCf^g* = 0. (139) 

The conservation of property (139) under time evolution allows us to introduce a 
notion of complex germ. 

Let otij be symmetric complex matrix n x n such that Ima > 0. Consider the 
following function g a G iS(R n ): 

g a (0 = exp ^ jr SkaM&j (140) 

The following definition is analogous to definition 2. 

Definition 7 [10]. The following n- dimensional subspace of the complex In- 
dimensional space R n : 

n 

Qa = {(Pl,-,Pn',Ql,-,Qn)\Pi = 

J'=l 

wra// 6e referred to as a complex germ corresponding to the matrix a. 
The following lemma is the analog of lemma 8. 
Lemma 20. 1. Let (p,q) G Q a . Then 

1 d , 
(PS - Q-^)9a = 

for arbitrary (p, q) G Q a . 

£.Le£ (p£ — q\j^)f = 0. TTien / = cg Q /or some constant c G C. 
The proof is straightforward. 

Let (5P , 5Q°) G £ Q o. Then (<5P*, G </ a t for the matrix a 1 being a solution to 
the Riccati equation: 

d t d 2 H d 2 H t t d 2 H t d 2 H t . . 

— a — a — a a a (141) 

dt dQdQ dQdP dPdQ OQdQ { ' 

Lemma 20 and eq.(139) tell us that the initial condition g° = g a o evolve into the 
following function 

9* = c l g a t. 

for some constant c* G C. Thus, we see that Gaussian wave packet (134) evolve into 
a Gaussian one as h — > 0. One can also consider complex germ creation operators 
of the form (139), where (<5P**, b~Q u ) G Ga* and find another asymptotic solutions 
to eq.(132). Thus, the presented approach allows us to find unambiguously complex 
germ creation and annihilation operators and the Riccati equation (141) without 
consideration of eq.(132). 
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Appendix B 



In this appendix we show how one can reproduce for the operator-valued case the 
classical Hamiltonian system (133), as well as the Hartree-like equation (66) by mak- 
ing use of the equations for mean values of the observables and of the argumentation 
analogous to the derivation of the Ehrenfest theorem (see, for example, [20]). Note 
also that the method to be presented in this appendix is heuristic. 

1. Consider the ordinary quantum mechanical Schrodinger equation 



in Sj^ = H ( Q ^ h s_^_ i s_ u , (Qiq) (142) 



corresponding to the physical system which is " semiclassical" with respect to Q G R" 
and "quantum" with respect to q G R m . Consider such solutions ip l £ L 2 (R n+m ) 
that 

(V, A (o, -ift A, Qi A - L\ A (q\ P\ q, A (143) 



for ^-independent functions A. For example, the mean values of "semiclassical" ob- 

'dQ> 



servables A(Q, —ih-jBj) are required to have limits as % — > 



Notice that there exist wave functions obeying the condition (143): examples of them 
are functions like 

where g f G S(R n+m ). 

We are to find possible equations for P*, Q 1 . To do this, consider first the equation 
for the mean value of some observable B = B(q, —i-j^): 

th-^ t ,B^) = ^ t ,[B,H]^ t ), (144) 

where H = H(Q, —ih-J^, q, As the left-hand side of eq.(144) is of order 0(h) as 

h — > 0, the right-hand side of it should also vanish in a leading order of h. Eq.(143) 
implies that one should replace Q by Q l and —ihd/dQ by P* for calculating the 
leading order of the right-hand side of eq.(144). 
Therefore, one should demand 

W,[B,HW) = 0, (145) 
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where H = H(Q\ P\ q, -ijfc). Eq.(145) implies that for arbitrary B Tr(B[H, II']) 

aplies th 



' dq 

0, where II* is a projector on ip t . This implies that [H, IT] = 0,i.e 



for some A = A(P*, Q f ) G R. 

Consider now the equation for the mean value of the observable A(Q,—ih-J^). 
Making use of eq.(143), one finds 

.dA ( . .( dA dH dA8H\ A _ 2 . 

where A = A(Q t ,P t ). We have taken into account that the semiclassical approxi- 
mation for the commutator of two "semiclassical" observables is [21] their Poisson 
bracket multiplied by ih. By using the relation 

(i/> t ,6Hi/> t )=6\ (146) 

one has 

dA OA d\ OA d\ 



dt dQtdPt dPtdQ 1 ' 

Therefore, we have reproduced the Hamiltonian system being used in constructing 
asymptotics for the operator- valued case [14]. 

2. An analogous technique can be also used for derivation of eq.(66) for the 
operator-valued case. Let us illustrate the approach for the case of constructing 
asymptotics for eq.(107). 

First of all, let us represent eq.(107) through the creation and annihilation opera- 
tors. Analogously to section 6, consider such Hilbert space Ti that R 2 ®L 2 (R ra ) C Ji. 
We choose Ti as a space of sets of functions 

$„, m : W x ... x R" x {1, 2} x ... x {1, 2} = K un x {1, 2} m -> C 

which satisfy the condition 

oo 2 . 

jdx 1 ...dx n \$ ntm (x 1 ,...,x n ,I 1 ,...,I m )\ 2 <oo 

n,m=0 7 1 ,...,7 m =l 

and are symmetric separately with respect to Xi and with respect to Jj. We identify 
such element $ G H that 

a) & n ,m — as n ^ N or m ^ 1; 

b) $„,i(xi, x N , I) = $ N ,i(xi, xn) 

with the element \I/ G R 2 ®L 2 (R ra ). By n we denote the projector on the subspace 
R 2 <g> L 2 (R ra ) of the Hilbert space H. 
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We introduce the following creation and annihilation operators in 7i: 

(a + (x)$) n:Tn (x 1: x n , h, I m ) = 

1 n 

/= y ' 8{x Xj)$ n _i m (Xi, Xi—i, X n , Ii, Im)-, 

V n i=i 

(a"(x)$) n _i >ro (a:i,...,a; n _i,/i,...,/ m ) = v 7 ^ h,..;Im), 

< ^)n,m(^l, X n , -Tl, An) = X! ^n.m-l (^1 j •••) ^n; A, -(j'+l) •••> An)j 

V m j=l 

(67$) 

where iGR",/6 {1,2}. 

Consider the operator H : H — > 7i of the form 

NH(a + /VN,a-/VN,b + ,b-) = dxa + (x)a~ (x) ]T biB u (x)bj + 

J i,j=i 

l r h 2 

-( / dxa + (x)[ A + C/"(x)]a"(a:)+ 

n. 7 2m 

1 /■ 2 
+7^j J dxdya + {x)a + {y)V{x,y)a (x)a (y))J2 b t b i- 

Analogously to lemma 10, one shows that any solution to eq.(107) obeys also the 
following equation 

d 

(i—-NH)¥ = 0. (147) 

Let us analyse this equation in a way analogous to the consideration of eq.(142). 
Consider such solutions to eq.(147) that 

A(a + /VN, cT/VN, b + , &")**)- 

- A(^\ <p\ b + , b~)¥) (148) 

for any polynomial function A being invariant under substitution a ± (x) — > a ± (x)e ±ia , 
bf — > bfe ±llS , a,(3 = const. Example of ^* obeying eq.(149) is 

^ = i^G- 

Let us derive now eq. (108). Consider the mean value of the operator B = B(b + , b~). 
One has: 

i^-i^^m 1 ) = N(¥, [B,H]¥). 
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As the right-hand side of this equation is of order O(N) as N — > oo, one should 
require it to vanish. Analogously to the quantum mechanical case, this implies that 

(JT(^y,6+ 6-)-A)«* = 0. 

Making use of the definition of the operator H, one shows that when ^ G R 2 <g) 
L 2 (R un ), there are two eigenvalues A, 

\ = H±(<p u ,<p t ,b + ,b-), 

and 

K,i,± = cH^*^ t )x t N ,±- 

These eigenvalues and eigenf unctions coincide with those which has been found in 
section 9. 

As the commutator between operators a ± /\/~N tends to zero as N — > oo, one can 
use ordinary semiclassical technique [21] to compute the commutator [A, if]. Let 



A = A{a + /VN,a'/VN). 

Then 

[A,H]¥)~ 

~ [ 'J \6<p(x) 5 V *(x) ~ 6<p(x) 6<p*(x)J ) ^ N ^°° ' 

we have omitted the arguments tp*, tp of the function A and the arguments <p*, ip, b + , b~ 
of the function H. Making use of eq.(146), one obtains 

^ = / & (J^_^_i^_^, (149) 

dt J \6ip(x)6ip*(x) 5<p(x) 5<p*(x) J 

where the sign + or — depends on the choice of the wave function We can notice 
that eq.(149) is a consequence of the Hartree-like equation (108). Thus, the approach 
based on the Ehrenfest theorem allows us to derive the equation for </?*. 
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